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1.0  INTRODUCTION 


Holography  is  well  known  for  its  fascinating  three-dimensional  images.  In  the 
scientific  community,  however,  holograms  are  frequently  used  for  obtaining  three- 
dimensional  data  such  as  size  and  shape  about  the  objects  of  which  a  hologram  is  made. 
A  common  technique  for  extracting  such  data  is  illustrated  in  Fig.  1.1.  This  approach 
involves  taking  pictures  of  the  optical  image  from  different  viewing  angles  and  storing 
them  in  a  computer.  Then,  the  computer  uses  the  stored  pictures  to  reconstruct  the  three- 
dimensional  image  of  the  object. 


Figure  1.1  Optical  reconstruction  of  holograms. 


Reconstruction  via  digital  holography  takes  a  different  approach.  The  steps 
involved  are  illustrated  in  Fig  1.2.  The  first  step  is  to  magnify  the  hologram.  Next,  an 
electronic  camera  and  a  frame  grabber  are  used  to  electronically  capture  and  digitize  the 
hologram  for  storage  in  a  computer.  Finally,  the  computer  is  used  to  mathematically 
reconstruct  the  image  of  the  object  from  the  fringe  pattern  in  the  hologram. 

The  impetus  for  exploring  digital  holography  came  as  a  natural  progression  from 
the  conventional  technique  depicted  in  Fig  1.1.  Since  a  computer  had  to  be  used  anyway, 
it  was  thought  that  if  the  hologram  itself  could  be  directly  digitized  into  a  computer,  then 
the  step  of  using  a  camera  to  take  pictures  from  various  angles  with  its  associated  error  of 
mechanical  placement  could  be  eliminated.  Additionally,  since  the  hologram  itself  would 
be  in  the  computer,  this  would  allow  processing  of  the  hologram  before  the  images  were 
formed.  The  recreation  would  be  a  matter  of  mathematics  rather  than  optics. 
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2.0  THEORETICAL  BACKGROUND 


Diffraction  is  the  key  theory  for  describing  how  electromagnetic  waves  or  light 
propagate  through  space.  Sommerfeld  described  diffraction  as  “any  deviation  of  light 
rays  from  rectilinear  paths  which  cannot  be  interpreted  as  reflection  or  refraction” 
[Goodman,  1996,  p.  33].  By  utilizing  the  principles  of  diffraction,  one  can  characterize  a 
wave  at  locations  along  its  path  of  propagation.  In  order  to  take  a  complete  snapshot  of  a 
wave  at  any  particular  location,  it  is  necessary  to  store  both  the  amplitude  and  phase  of 
the  wave.  A  recording  of  this  type  is  known  as  a  hologram  which  means  a  “total 
recording”  [Goodman,  1996,  p.  295]. 

2.1  Mathematics  of  Holography 

The  purpose  for  presenting  the  mathematical  theory  of  holography  is  to  allow  the 
quantifying  of  the  holographic  process  so  that  its  limitations  can  be  predicted. 
Development  begins  with  an  object  beam  and  illustrates  how  approximations  to  the 
diffraction  theory  can  be  applied  to  simplify  description  of  the  object  beam  as  it 
propagates  through  free  space.  Eventually,  the  approximations  will  enable  the 
description  of  the  object  beam,  after  propagating  through  some  distance  d,  in  terms  of  the 
Fourier  transform.  Next,  the  method  of  recording  the  object  beam,  both  amplitude  and 
phase,  will  be  discussed.  Finally,  the  process  of  optical  reconstruction  of  the  object  beam 
will  be  explained  thus  resulting  in  the  recovery  of  the  original  object  waveform.  By 
laying  this  theoretical  background  at  this  stage,  we  can  then  use  this  background  to 
predict  the  benefits  and  limitations  of  reconstructing  the  beam  digitally. 

2.1.1  The  Object  Beam 

Consider  Fig.  2.1  which  contains  an  infinite  opaque  screen  with  an  aperture.  It  is 
through  this  aperture  that  the  object  wavefront  or  beam  originates.  Additionally,  a  region 
of  observation  is  shown.  The  region  of  observation  is  assumed  to  be  a  plane  parallel  to 
the  original  screen.  Beginning  with  the  Rayleigh-Sommerfeld  diffraction  formula,  the 
object  beam,  in  the  observation  plane,  can  be  described  by  the  Huygens-Fresnel 
relationship  [Goodman,  1996,  p.  46-66] 

U(£,  v)  =  JJ  g(£>  n;  *,  y)0(x,  y)dxdy ,  (2.1) 

I 

where  0(x,  y)  is  the  object  beam  at  the  screen  and  g(^,r];x,  y)  is  a  transformation  kernel 
[Yaroslavskii,  1980,  p.  4].  In  Eq.  (2.1),  the  transformation  kernel  is  given  by 

1  ejkrPP' 

g(Z,mx,y)  =  — - cos(0).  (2.2) 

A  rpp' 

0(x,y)  is  zero  outside  the  aperture  X-  In  Fig.  2.1,  P'  represents  a  point  in  the 
observation  plane  for  which  the  wavefront  is  being  calculated  and  P  represents  one  point 
in  the  aperture  over  which  the  integration  is  being  carried  out.  The  angle  6  is  that 
between  the  line  connecting  the  points  P  and  P'  and  the  normal  of  the  aperture  plane. 
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Figure  2.1  Diffraction  geometry. 


Letting  rpp>  represent  the  line  between  P  and  P'  and  d  the  distance  between  the  planes, 
then  from  the  geometry, 

cos(0)=  — .  (2.3) 

rPP' 

Substituting  Eq.  (2.3)  into  Eq.  (2.2)  yields 

d  eik'PP' 

3(4, n  \x,y)  =  — - 5-.  (2.4) 


Approximations  can  be  made  to  reduce  the  formula  into  a  more  usable  form 

2 

[Goodman,  1996,  p.  66-67].  If  0  is  small,  then  (rpp>)  in  the  denominator  can  be 

2 

approximated  as  d  .  However,  this  approximation  cannot  be  made  in  the  exponential 

term  of  the  numerator  since  small  changes  in  rpp>  will  cause  large  variations  in  e  pp  . 

Approximations  can  be  made  in  that  term  through  use  of  the  binomial  expansion  of  a 
square  root  as  follows.  The  distance  between  P  and  P'  is  given  exactly  by 


rpp,='yl^  +(£-*)  +(v~y)  -dJi  + 


(£-*) 


v-y 


The  binomial  expansion  of  a  square  root  is 

Vl+F  =  1  +  -b--b2+...  (1 

2  8 

for  1^1  <  1 .  Assuming  the  1st  two  terms  will  give  an  adequate  approximation,  we  have 


rpp'Bdl+i{ij£) +i[27z)  ■ 

Substituting  this  approximation  into  Eq.  (2.4)  and  rearranging  terms  produces 

(  jkd  \  j* 

9  (%,ri;x,y)=  -j^-je2d^-x)2  +  (ri-y)2^ 

which  leads  to 

U(£,t7) = Jj  °(*>  ^exP  +(*?-  y)2) 


=  -r^-exp  ^(£2+772)  JJO(x,y)exp  ^-( *2  +  y 2)  exp  - 
7  L  Jj-  L  J  L 


■(x^  +  yrj)  uxdy. 


Equation  (2.9)  is  the  Fresnel  approximation  relating  the  complex  optical 
amplitude  distribution  in  the  observation  plane  to  that  at  the  screen  [Goodman,  1996,  p. 
67].  When  d  is  large  enough  for  the  approximation,  the  observer  is  in  the  Fresnel  or  near¬ 
field  diffraction  region.  In  the  last  exponential  term,  the  wave  number  k  was  substituted 
with  its  equivalent  2 n  I A .  Notice  that  the  integral  in  Eq.  (2.9)  is  just  the  two- 

jk(x2  +  y2) 

dimensional  forward  Fourier  transform  of  0(x,  y)  exp -  ,  except  for  the 

2d 


factor  —  in  the  last  exponential  term.  We  will  explore  this  factor  later  when 
Ad 

investigating  digital  reconstruction. 

The  significance  of  this  development  can  begin  to  be  appreciated.  Computing  a 
Fourier  transform  digitally  on  a  computer  can  easily  be  done  through  use  of  the  fast 

T  / k(x2  +  y2)l 

Fourier  transform  (FFT).  Also,  the  term  exp  -  in  Eq.  (2.9)  can  be  created 

2d 

mathematically  as  can  all  of  the  terms  in  front  of  the  integrals.  So,  if  the  wavefront 
0(x,  y)  in  Eq.  (2.9)  is  known,  then  can  be  computed  mathematically  on  a 

computer. 

A  further  simplification  to  Eq.  (2.9)  can  be  developed  and  is  known  as  the 
Fraunhofer  or  far-field  diffraction  approximation.  If  d  is  very  large  such  that 

J»k(x2  +  y2)max 


jk(x2+y 2) 


(2.10) 


which  leads  to 
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ejkd 

u(^,T?)=7a7exPi 


#«2+"2) 


jj°^  ,j)exp 


J2k_ 

M 


(x<i;  +  ytj)  pcrfy  (2.11) 


Thus,  in  the  Fraunhofer  region,  the  optical  complex  amplitude  is  related  to  0(x,  y)  by 

the  forward  Fourier  transform  directly  (except  for  the  factor  —  in  the  last  exponential 

Ad 

term  as  mentioned  earlier). 


2.1.2  Recording  a  Complex  Amplitude  Hologram 

To  record  a  hologram,  it  is  necessary  to  store  information  about  both  the 
amplitude  and  phase  of  the  object  wavefront.  The  standard  technique  to  accomplish  this 
is  by  recording  the  constructive  and  destructive  interference  between  two  coherent 
wavefronts;  this  process  is  known  as  interferometry.  The  mathematical  description  of 
this  technique  can  be  illustrated  by  computing  the  superposition  of  two  wavefronts.  The 
wavefront  from  the  object  beam,  at  the  film  plane,  can  be  described  as 

U(^,77)  =  C/(^,77)e-^’T?),  (2.12) 

where  and  0(£,7 7)  are  the  amplitude  and  phase  of  the  object  beam,  respectively. 

Similarly,  the  reference  beam  is  given  by 

=  (2.13) 

where  A  and  y/(^,rj)  are  the  amplitude  and  phase  of  the  reference  beam,  respectively. 

By  design,  the  reference  beam  amplitude  is  made  approximately  a  constant. 

Unfortunately,  photographic  film  cannot  record  complex  values.  Film  actually 
records  intensity,  which  is  calculated  by  multiplying  one  term  by  its  conjugate.  For  the 
case  where  there  are  two  wavefronts  incident  on  the  film  simultaneously,  the  film  will 
record  the  intensity  of  the  sum  [Kozma,  1966],  Thus,  the  intensity  of  the  object  beam  and 
reference  beam  together  is 

/(£,77)  =  (U+A)(U  +  A)* 

=  |U|2  +|A|2  +  A*U+ AU*  (2.14) 

=  A2  +  U  2  (£,  77)  +  2  A  (/  (£,  77)  cos[V  (£,  77)  -  0(f ,  77)] 

The  last  term  of  Eq.  (2.14)  shows  how  both  amplitude  and  phase  of  the  object  are 
preserved.  In  contrast,  if  the  wavefronts  incident  upon  the  film  are  not  coherent,  then 
constructive  and  destructive  interference  will  not  occur.  In  this  case,  the  film  will  record 
purely  the  addition  of  the  individual  intensities  as  in  ordinary  photographs  [Kreis,  1996, 

p.  22], 


Film  does  not  necessarily  record  intensity  linearly.  However,  film  does  typically 
have  a  linear  portion  in  its  t-E  curve.  If  the  exposure  time  with  respect  to  intensity  is 
chosen  so  that  exposure  is  in  the  middle  of  this  linear  portion  and  if  A  is  large  enough 
compared  to  U (E,,vj) » then  the  recording  will  be  confined  to  the  linear  portion  of  the  t-E 
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curve  [Kozma,  1966].  Assuming  this  linearity  in  the  mapping  of  intensity  exposure  to 

r\ 

amplitude  transmittance  and  that  |A|  is  uniform  across  the  recording  film  surface,  then 
the  transmittance  is 

r(£,fj)  =  h  +P'(\V\2  +  A*U  +  AU*),  (2.15) 

where  tj,  is  a  uniform  bias  resulting  from  the  reference,  and  P'  is  the  product  of  the 
slope  (3  of  the  film’s  t-E  curve  at  the  bias  point  and  the  exposure  time  [Goodman,  1996, 
p.  303].  /?'  is  negative  for  a  negative  transparency  and  positive  for  a  positive 
transparency. 


2.1.3  Image  Reconstruction 

The  development  in  this  section  will  follow  closely  with  [Goodman,  1996,  p.  302- 
309].  Illuminating  the  hologram  by  a  coherent  reconstruction  wave  B(£,7])  produces 
light  transmitted  by  the  transparency  given  by: 

B(f ,  77>r(&  77)  =  t*B  +  j0OJU*B  +  P' A*BU + /TABU*  =  Sj  +  S2  +  S3  +  S4.  (2.16) 

If  B  =  A ,  then  S3  =  /T  A  U  which  is  a  constant  times  a  duplicate  of  the  original  object. 

If  B  =  A* ,  then  S4=/3'|A|2U*  which  is  a  constant  times  the  conjugate  of  the  original 

object.  Optically,  if  the  reference  wave  used  in  the  reconstruction  process  is  the  same  as 
the  original  reference  wavefront,  then  a  duplicate  of  the  original  diverging  object 
wavefront  will  be  generated.  This  will  produce  a  virtual  image  that  can  be  imaged  onto  a 
camera  through  use  of  a  lens.  On  the  other  hand,  if  the  conjugate  of  the  original 
reference  wavefront  is  used,  then  a  real  image  is  created  corresponding  to  an  actual 
focusing  of  light  and  so  use  of  a  lens  is  not  necessary  [Goodman,  1996,  p.  301]. 


Dennis  Gabor  is  accredited  with  the  first  appearance  of  holography  in  the  1940s. 
However,  it  was  work  with  the  off-axis  reference  wave  by  Leith  and  Upatnieks  in  the 
early  1960s  and  the  invention  of  the  laser  that  really  led  to  widespread  research  in  the 
area  of  holography  [Reynolds,  1989,  p.293].  Holograms  using  an  off-axis  reference 
wave  are  generally  referred  to  as  Leith-Upatnieks  or  offset-reference  holograms.  It  is 
illustrated  in  Fig.  2.2.  It  is  made  by  introducing  a  separate  distinct  reference  wave  at  an 
offset  angle  to  the  object-film  axis.  For  this  case,  the  reference  wave  is  described  by 


A({,17)  =  Ae-^ . 


(2.17) 


where  the  spatial  frequency  a  of  the  reference  wave  along  the  film  plane  is  given  by 

a  =  ^-.  (2.18) 


Thus,  the  intensity  distribution  across  the  film  for  the  Leith-Upatnieks  hologram  is 
/(|, 77)  =  (Ae-i2™*  +U&T1)e-j^){Ae-i2nar' 

=  A2  + 1/  2  (£,  77)  +  2  A(7(£,  77)  cos(2jtar]  -  </>(?;,  77)) 


(2.19) 


Figure  2.3  is  an  actual  scanned  hologram.  It  serves  to  illustrate  what  is  actually 
recorded  on  by  a  hologram.  The  line  structure  that  can  be  seen  is  usually  referred  to  as 


« 


m 
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the  fringe  pattern.  Note  that  the  picture  is  an  enlargement.  Ordinarily,  the  fringe  pattern 
on  a  hologram  cannot  be  seen  with  the  eye  unless  viewed  through  a  microscope. 


M 


film 


Figure  2.2  Leith-Upatnieks  hologram  setup. 


Figure  2.3  Picture  of  an  actual  hologram  as  seen  through  a  microscope. 


After  the  film  is  developed,  its  resulting  transmittance  is 

T&ri)  =  tb  +  j3'(|U|2  +  AUejl7caT1  +  AU*e-j2KCC11). 
This  can  be  reduced  to  the  following  four  terms: 


(2.20) 
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*1  =<b 


h  =/J'|U(^,r,)|2 

(2.21) 

t3  =  P'AV(Z,ri)ej2na11 
\A=p'AU\^i])e~j2naj] 

Figure  2.4  illustrates  the  optical  reconstruction  of  such  a  hologram.  Illumination  by  a 
normally  incident,  uniform  plane  wave  of  amplitude  B,  which  is  a  duplicate  of  the 
original  reference  beam,  gives 

s  i=tbB 
s2  =p'B\\Mtnf 

(2.22) 

53  =  p'BA\K^li)ei2naT] 

54  =  p'BAV*(%,Ti)e~j2nan 

Sj  is  an  attenuated  version  of  the  incident  reconstruction  field  and  so  represents  a  plane 
wave  traveling  down  the  transparency  axis.  S2  is  spatially  varying  and  so  contains  plane 
wave  components  traveling  at  various  angles.  S3  is  proportional  to  the  original  object 
wavefront  U  which  implies  the  term  generates  a  virtual  image  at  distance  d  and  angle  0 
from  the  transparency.  S4  is  similar  to  S3 .  It  is  proportional  to  the  conjugate  of  the 


Figure  2.4  Images  reconstructed  from  a  Leith-Upatnieks  hologram  [from  Goodman, 
1996,  p.  307], 


9 


object  wavefront  U*  and  produces  a  real  image  on  the  opposite  side  of  the  virtual  image 

at  an  angle  -0  to  the  transparency  axis.  These  results  indicate  that  both  the  real  and 
virtual  images  are  angularly  separated  from  each  other  and  the  waves  S3  and  S4  .  This 
is  due  to  the  angular  offset  of  the  reference  wave.  Also  indicated  is  that  an  exact 
duplicate  or  conjugate  of  the  original  reference  wave  is  not  always  needed  for 
reconstruction. 

2.1.4  Separation  of  Components  in  Reconstruction 

In  order  to  obtain  separation  of  the  twin  images  and  the  light  transmitted  along  the 
transparency  axis,  the  spatial-frequency  content  of  Sj  through  S4  above  should  not 
overlap.  This  is  achieved  through  appropriate  selection  of  the  reference-offset  angle  9 . 
Determination  of  the  minimum  angle  can  be  found  from  the  Fourier  transform  of  Sj 
through  S4  ,  which  will  be  denoted  as  Fj  through  F4 .  The  transform  of  the  object 
wavefront  will  be  denoted  as  Fu .  We  will  assume  the  highest  frequency  content  of  the 
object  is  W.  From  these  spectra,  the  spatial  carrier  frequency  a  needed  to  achieve  image 
separation  is  found. 

Figure  2.5  illustrates  these  various  spectra.  The  spectrum  Fj  is  an  impulse  at 
origin  and  thus  has  zero  bandwidth.  The  spectrum  F2  is  proportional  to  the 
autocorrelation  of  Fu  and  so  its  bandwidth  is  2 W.  The  spectrum  F3  is  proportional  to  Fu 
but  displaced  to  center  frequency  a  and  so  its  bandwidth  is  W.  The  spectrum  F4  is 
proportional  to  a  reflected  version  of  Fu  and  is  centered  at  -a.  Its  bandwidth  is  also  W. 
From  the  graphs,  |F3|  and  |F4|  can  be  isolated  from  |F2|  if  a  >  3W ,  or  sin#  >  3 WX  . 

Thus, 

0min  =  sin-1  3 WX  .  (2.23) 

If  the  reference  wave  is  much  stronger  than  the  object,  then  U  «  A  and  so  |F2|  is  much 
smaller  than  |F[|,  |F3|,and  |F4|.  In  this  case. 


Figure  2.5  Spectra  of  components  from  a  reconstructed  Leith-Upatnieks  hologram. 
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^min  =  sin-1  WX .  (2.24) 

Note  that  Fig.  2.5  indicates  frequency  displacement  in  the  y-dimension  only  since  the 
reference  wave  was  offset  only  with  respect  to  the  y-axis.  Similar  analysis  would  apply 
for  an  offset  angle  with  respect  to  the  x-axis. 

2.2  Spatial  Frequency  Due  to  Interference 

The  issue  addressed  here  is  to  develop  a  formula  for  the  spatial  frequency  that  is 
recorded  on  the  photographic  film  used  to  record  the  hologram.  This  is  necessary  for 
determining  the  highest  spatial  frequency  recorded  on  the  hologram,  which  will  then 
dictate  the  Nyquist  sampling  frequency  required  to  avoid  aliasing.  The  spatial  frequency 
results  from  the  coherent  interference  between  the  reference  wave  and  the  object  wave. 
Figure  2.6  illustrates  the  geometry  useful  for  analyzing  this  phenomenon.  As  will 
become  apparent  later  in  the  development,  it  is  the  maximum  angle  between  the  object 
and  the  reference  wave  that  determines  the  highest  spatial  frequency  that  will  be  recorded 
for  any  particular  setup.  The  angles  indicated  in  Fig.  2.6  represent  the  maximum  for  the 
geometry  in  the  figure. 


object 


reference 


Figure  2.6  Object  and  reference  wave  angles. 

As  already  discussed,  the  film  will  record  the  intensity  of  the  superposition  of  the 
object  and  reference  waves  at  each  point  on  the  film.  The  intensity  is  given  by 
multiplying  the  sum  of  the  two  waves  by  the  conjugate  of  their  sum.  For  this 
development,  the  most  general  form  of  representing  waves  is  used  to  show  the  full 
mathematical  development. 
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Defining  the  intensity  as  /(^,  77) ,  the  object  wave  at  the  hologram  plane  as 
U(£,77) ,  and  the  reference  wave  as  A(^, 77) ,  then  I(%,ri)  is  given  by  the  following: 

I(Z,ri)  =  [  A(£,  77)  +  U(£,  7])]  •  [Atf,  77)  +  U(|,  77)]* 

=  [  A  exp-  ;'(4  •  r  -  wt  +  (pa  )  ■ +  U  exp-  j[ku  ■  r  -  wt  +  <j)u )]  (2.25) 

|  A  exp-  j(ka-r -wt  +  <t>a)  +  U  exp-  j(ku  •  r  -  wt  +  <j)u  )]  . 

This  simplifies  to 

/ (<lj, 77)  =  2Acos(ka  -  r  —  wt  +  <payu  cos(ku  •  r  -  wt  +  +  A2 


From  Fig.  2.6, 
and 


+2Asin(fca  -  r  -wt  +  <j)ayu  sin(&H  -  r  -wt  +  (j)u^  +  U^ 

=  2AU cos[(&a  r  -  wt  +  ^a^-{ku  r  -  wt  +  </>M)j  +  A2  +C/2. 

A  A 

=  E,k  sin(0a )  +  Qi  cos  {9a ) 


(2.26) 


(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


,  kU  =  Sin(^M  )  +  Ck  COS  (6U)  . 

Since  the  position  vector  is  given  by 

r  =U  +  rjri  +  &, 

then 

ka  r  =%k  sin(0a )  +  £k  cos (9a ) 
and 

ku  -r  =  -£k sin(0 u)  +  Qc cos(0 u) . 

Substituting  Eq.  (2.30)  and  Eq.  (2.31)  into  Eq.  (2.26)  and  using  £  =  0  for  the  hologram 
plane  leads  to 

/(£,  77)  =  2AU cos{[/:£ sin(0fl  )  -  wt  +  <f>a ]  -  [-&£ sin(0K )  -  wt  +  ]}  +  A2  +  U2 

(2.32) 

=  2AUcos{fc£[sin(0a)  +  sin(0u)]  +  0a  -0„}  +  A2  +U2. 

Since  (j)a  and  (pu  are  constant  phase  terms,  then  the  spatial  frequency  a  of  the  resulting 
fringe  pattern  can  be  determined  from 

2nat;  =  k^\sm{9a )  +  sin(0„ )].  (2.33) 

2  jt 

Using  k  =  —  and  solving  for  a  in  Eq.  (2.33)  gives  the  resulting  formula 

A 

sin(0a)  +  sin(0„) 

2 

Note  that  for  a  given  total  angle  between  the  reference  and  object  waves,  the 
spatial  frequency  a  can  still  vary.  For  example,  let  6a  =  30°  and  9U  -  60° .  Then 


(2.34) 
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a  =  However,  if  0a  =  45°  and  0M  =  45° ,  then  a  =  — - — .  So,  even  though 

X  X 

0a  +  eu=  90°  in  both  cases,  a  is  different  for  each  case.  It  can  be  shown  that,  for  a 

given  sum  0a+9u,  the  maximum  a  occurs  when  9a=6u.  Figure  2.7  illustrates  the 

relation  between  a,  X,  6 a,  and  0U. 


Figure  2.7  Plot  of  spatial  frequencies  resulting  from  interference  between 
various  reference  plane  wave  angles  0a  and  object  angles  0U  for  X=632.8nm. 
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3.0  DIGITAL  RECONSTRUCTION 


Digital  reconstruction  of  holograms  has  several  advantages.  One  is  that  it  avoids 
the  speckle  that  occurs  during  the  optical  reconstruction  (but  not  the  speckle  recorded  on 
the  hologram).  Another  advantage  is  that  the  raw  data  recorded  on  the  hologram  is 
utilized  and  so  allows  more  direct  reconstruction  of  the  image.  Another  advantage  is  that 
the  depth  of  field  is  large  which  is  good  for  reconstruction  as  in  tomographic  back- 
projection. 

Of  course,  digital  reconstruction  has  its  disadvantages  as  well.  Some  of  its 
limitations  are  addressed  in  this  chapter.  For  instance,  the  large  depths  of  field  are  a 
result  of  low  resolution.  Each  of  the  steps  shown  in  Fig.  1 .2  has  an  effect  on  the  quality 
of  the  reconstructed  image. 

3.1  Digitization  of  the  Hologram 

3.1.1  Microscope  Objectives  and  Numerical  Aperture 

In  order  to  scan  a  hologram  into  a  computer,  the  instrument  must  resolve  the  very 
small  fringes  the  are  recorded  on  the  hologram.  Typically,  it  is  necessary  to  utilize  a 
microscope.  So,  the  attention  is  turned  to  the  ability  of  microscopes,  or  more 
appropriately,  the  microscope  objective  to  resolve  the  fringes. 

The  mathematical  development  starts  with  the  Fraunhofer  approximation  for  the 
field  amplitude  given  in  Eq.  (2.1 1).  The  integral  in  that  equation  is  carried  out  over  the 
aperture  of  the  system.  For  the  case  of  a  microscope  objective,  the  aperture  is  circular. 
Using  a  circular  aperture  as  the  limits  of  integration  and  assuming  a  plane  wave  incident 
upon  this  aperture,  Eq.  (2.1 1)  develops  into  a  Fourier-Bessel  transform.  The  solution  of 
this  transform  equation  is  well  known  and  the  resulting  intensity  described  is  known  as 
the  Airy  pattern  [Hecht,  1998,  p.459-461]. 

The  Airy  pattern  is  the  physical  phenomenon  that  is  used  to  establish  the 
resolution  limit  of  an  optical  system.  When  this  limitation  applies,  the  optical  system  is 
said  to  be  diffraction-limited.  Depending  on  the  quality  of  the  optical  components  in  the 
system,  the  resolving  power  may  be  considerably  worse  than  the  idealized  diffraction- 
limited  case. 

Assuming  a  diffraction-limited  optical  system,  a  quantifiable  measurement  of 
the  resolving  power  of  a  system  is  needed.  One  method  used  is  the  Rayleigh  criterion. 
This  criterion  states  that  two  points  can  be  resolved  if  the  maximum  of  one  point’s  Airy 
pattern  overlaps  at  the  first  zero  of  the  other  point’s  Airy  pattern  [Longhurst,  1967,  p. 
302].  This  is  illustrated  in  Fig  3.1.  From  the  superposition  principle,  the  Airy  pattern  of 
each  point  adds  to  give  the  result.  In  Fig  3.1,  there  is  a  line  that  droops  down  between  the 
peaks  of  the  individual  Airy  patterns.  That  droop  indicates  that  the  individual  points  of 
light  can  still  be  distinguished. 
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Summation  of  both 
Airy  patterns 


Figure  3.1  Illustration  of  the  Rayleigh  criterion. 


Using  the  Rayleigh  criterion,  the  limit  of  resolution  is  then  the  radius  to  the  first 
dark  ring.  This  radius  is  given  by 


,,  0.6  ir 

n  =- 


sin  6'  ’ 

where  X  is  the  wavelength  in  the  image  space  and  6'  is  the  semi-angle  of  the  emergent 
cone  of  light  [Longhurst,  1967,  p.  220-223].  Since 


(3.1) 


x'=h- 


n 


(3.2) 


where  A0is  the  wavelength  in  a  vacuum  and  n  is  the  refractive  index  of  the  image  space 
[Hecht,  1998,  p.  103],  then  Eq.  (3.1)  can  be  written  as 

0-61A, 


h' 


n'  sin#' 


(3.3) 


[Longhurst,  1967,  p.  224].  It  should  be  noted  that  Longhurst  [1967,  p.  224]  used  X 
instead  of  Xa  to  represent  the  wavelength  in  a  vacuum. 

The  Rayleigh  criterion  is  perhaps  the  most  commonly  used  criterion.  However,  at 
least  one  other  deserves  mention.  It  is  called  the  Sparrow  criterion.  The  Sparrow 
criterion  defines  the  resolution  limit  as  occurring  when  the  intensity  sums  result  in  a  flat 
top  or  plateau  between  the  two  peaks  versus  the  droop  shown  in  Fig  3.1  [Reynolds,  et.  al., 
1989,  p.  40-41],  It  can  be  inferred  from  Reynolds,  et.  al.  [1989,  p.  42]  that  the  Sparrow 
criterion  gives  a  resolution  limit  of 


y  _  0-473 6Ap 
n'  sin0' 


(3.4) 


When  using  a  microscope  (as  in  our  situation  of  scanning  the  fringes  from 
holographic  film),  the  resolving  power  requires  additional  development.  Figure  3.2 
depicts  the  setup  for  a  microscope.  In  this  case,  the  optical  sine  theorem  should  be 
applied  to  relate  the  image  points  to  the  object  points.  Using  Fig.  3.2,  this  relation  is 
given  by 


nh  sin  Q  =  n'h'  sin  9' , 


(3.5) 


r 


where  n  and  n  are  the  indices  of  refraction  of  the  object  space  and  image  space, 
respectively  [Hecht,  1998,  p.  264],  [Longhurst,  1967,  p.  297-298].  Equations  (3.3)  and 
(3.5)  can  be  combined  to  give  the  distance  h  between  the  object  points  at  the  diffraction 
limit  to  be 


h=0OK 

wsin# 


(3.6) 


[Longhurst,  1967,  p.  303].  The  quantity  nsin0  is  defined  as  the  numerical  aperture  and 
denoted  by  NA  [Longhurst,  1967,  p.  303].  Substituting  this  into  Eq.  (3.6) 
gives  the  resolution  limit  in  the  object  plane  for  a  diffraction-limited  system  as 


h=omK 

NA 


(3.7) 


Figure  3.2  Illustration  for  the  Optical  Sine  Theorem. 

3.1.2  Finite  Pixel  Size 

Ideally,  sampling  of  the  hologram  for  digitization  is  done  with  perfect  impulses. 
However,  this  cannot  be  done  with  real  components.  The  finite  width  of  the  true 
sampling  function  will  affect  the  reconstructed  image.  This  section  will  explore  the 
effect. 


After  the  hologram  is  magnified  with  a  microscope  of  appropriate  resolving 
power,  the  next  step  is  to  use  some  device  to  convert  the  optical  image  of  the  hologram 
into  an  electronic  form  that  can  be  stored  and  processed  in  a  computer.  A  typical  device 
used  is  an  electronic  camera.  The  component  in  the  camera  that  does  the  actual 
conversion  from  photonic  to  electronic  energy  is  referred  to  as  an  imager.  The  most 
common  imager  used  is  the  charge-coupled  device  or  CCD.  The  CCD  consists  of  an 
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array  of  pixels,  each  of  which  provides  a  packet  of  charge  to  the  camera,  which  converts 
the  charge  into  a  voltage  signal.  The  amount  of  charge  from  each  pixel  depends  on  a 
number  of  factors  including  the  amount  of  photon  energy  it  receives,  the  time  of 
exposure,  and  the  size  of  the  pixel.  It  is  the  size  of  the  pixels,  the  spacing  between  pixels, 
and  the  effect  they  have  on  the  image  quality  that  now  will  be  addressed. 


The  charge  output  from  each  pixel  is  the  result  of  accumulating  photon  energy 
over  its  area  and  converting  the  sum  of  this  energy  into  one  charge  packet.  Thus,  the 
output  is  actually  an  integration  of  the  energy  over  the  pixel  area.  We  make  the 
assumption  that  the  pixels  are  square  and  denote  the  length  of  each  side  as  W.  Also,  the 
pixel-to-pixel  spacing  in  the  £  and  77  dimensions  are  defined  as  %s  and  T)s ,  respectively. 
Figure  3.3  illustrates  these  dimensions.  Then  in  mathematical  terms,  the  charge  from  the 
mnth  pixel  is 


S.,  =  ]]u(ln)rect\ 


w 


lrect{} TT^77- 


(3.8) 


where  the  function  rect(  )  is  defined  by  [Kreis,  1996,  p.  268] 

,  ..  1 1,  \x\  <  — 

rect(x)  =  \  1  2 

0,  elsewhere. 


(3-9) 


When  the  charge  from  all  pixels  in  the  CCD  are  put  together,  the  result  is  a  sampling  of 
the  image  described  by 


m  n  —00—00  V 


£zm!L 

W 


rj-nr] 

W 


(3.10) 


■S(£-m%s)S{rl-nrls)- 


V 


Figure  3.5  Pixel  Size  and  Spacing. 
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It  will  now  be  shown  that  Eq.  (3.10)  is  equivalent  to  the  ideal  sampling  of  the 
convolution  of  the  image  with  a  rect  function  of  dimensions  the  same  as  a  pixel.  The 
continuous  convolution  [Ziemer,  et.  al.,  1983,  p.  44]  of  the  entire  image  with  a  2D  rect 
function  with  a  width  of  W  in  each  direction  is 


U.i&'n)  =  1  ]'U(a,b)rect(-^\ect 


r\-b 

~W 


dadb , 


(3.11) 


Then,  if  the  image  described  by  Eq.(3.1 1)  is  perfectly  sampled,  the  result  can  be 
described  by 


^f\-eJX~yadb5(i  -mDSfr-nn,).  (3.12) 


The  impulse  functions  in  Eq.  (3.12)  are  not  a  function  of  the  integration  variable  and  can 
be  brought  inside  the  integral.  This  gives 


Ucs(Z’Tl)  =  X£  J  \U{a,b)rect\  ^ -  \rect\  ]8{£-m%s)8{r\-nr\s)dadb .  (3.13) 


rj-b 


W 


W 


Using  the  fact  that  [Ziemer,  et.  al.,  1983,  p.  21] 

x(t)5(t-t0)  =  x(t0)8(t-t0) 

results  in 


ua(&v)  =  i  U  (a,b)rect 

m  n  “ 


\rect 


nris-b\ 
W  J 


dadb 


(3.14) 


(3.15) 


The  double  integral  of  Eq.  (3.15)  can  be  recognized  as  a  convolution.  Since  rect 
is  an  even  function,  Eq.  (3.15)  can  also  be  written  as 

'b-nr]s 


=  2,^1  \U(a,b)rect 


rect 


W 


\dadb 


(3.16) 


-8{Z-mZ,)8{r\-nr\t). 

Comparison  of  Eqs.  (3.10)  and  (3.16)  shows  that  they  are  equal.  Therefore,  it  can  be 
concluded  that  the  finite  pixel  width  and  its  operation  as  an  integrator  has  the  same  effect 
as  convolution  of  the  entire  image  continuously  with  a  rect  function  the  size  of  a  pixel 
followed  by  perfect  sampling. 


One  property  of  the  Fourier  transform  is  that  convolution  in  one  domain  is  the 
equivalent  of  multiplication  in  the  other  domain  [Ziemer,  et.  al.,  1983,  p.  132].  Therefore 


U(a,b)* 


rect 


(zL 


■a 


V  W 


3[t/(a,&)]3 


rect\ 


( mL 


■a 


y  w 


\rect\ 


nrh  zK 

W  . 


(3.17) 


where  *  indicates  convolution,  indicates  a  Fourier  transform  pair,  and  3  indicates 
the  Fourier  transform.  The  Fourier  transform  of  the  rect  function  is 


rect 


Ms 


w 


=  Wsinc(fW) . 


(3.18) 
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Assuming  the  spectral  extent  of  U(a,b)  to  be  ±FU ,  then  the  Fourier  transform  of  Eq. 
(3.15)  is  as  illustrated  in  Fig.  3.4.  Figure  3.4b  illustrates  the  special  case  where  the  pixel 

/ 

width  W  and  the  pixel  spacing  are  equal.  For  this  case,  the  amplitude  reduction  at  ~r 


due  to  the  sine  function  is  given  by 


WsincfW  = 


WsmjjtfW) 

7tfW 


=h. 

2 


2ff 

(*X^Xi)  * 


(3.19) 


This  corresponds  to  a  1.96  dB  reduction  in  frequency  amplitude  for  frequencies  at  half 
the  sampling  frequency. 


v 


Figure  3.4  Frequency  spectrum  of  image  convolved  with  a  pixel  of 
finite  width  and  perfectly  sampled:  (a)  General  case,  (b)  Special 
case  where  the  pixel  width  and  pixel  spacing  are  equal. 


3.2  Effect  of  Limited  Scan  Area 

The  width  of  the  pixels  is  not  the  only  effect  imposed  by  an  electronic  camera. 
Another  consideration  is  the  number  of  pixels  in  each  dimension.  The  digitized  hologram 
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obtained  from  a  camera  will  be  comprised  of  this  number  of  pixels.  As  a  result,  only  a 
small  portion  of  the  whole  hologram  is  digitized.  The  effect  of  this  is  examined  below. 

3.2.1  Maximum  Theoretical  Resolution 


Figure  3.5  illustrates  the  coordinate  system  for  the  object,  hologram,  and  image 
planes.  After  some  mathematical  development,  the  wavefront  at  the  observation  region 
(the  hologram  plane)  was  determined  to  be 


U(fr?)  =  — rexp 


jk_ 

2d 


jAd 
JJO(*,y)exp 


Of  +v2) 


] 

^(*2  +  /)Jexp|^-^(x£  +  y77)  dxdy. 


(3.20) 


Figure  3.5  Coordinate  system  for  the  object,  hologram,  and  image  planes. 
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The  goal  is  to  recreate  the  wavefront  from  the  object,  0(x,  y) .  The  integral  is  just 
the  two  dimensional  forward  Fourier  transform  of 


0(x,y)exp 


jk(x2 +  y2) 


2d 


(3.21) 


except  that  the  Ad  in  the  last  term  of  (3.20)  is  not  in  the  transform  definition.  This  must 
be  taken  into  account  for  proper  scaling  in  the  observation  plane.  This  can  be  accounted 
for  through  use  of  the  scale  change  theorem  of  Fourier  transforms.  The  scale  changes 
made  are 


P  = 


± 

Ad 


(3.22) 


and 


Substituting  these  into  Eq.  (3.20)  gives 


(3.23) 
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U(pAd,qAd)  = - exp 

jAd 


£((pAd)2+(qAd)2) 
2  a 


Jj  0(x,  y)  exp  (x2  +  y2 )  exp[-j2x(xp  +  yq)}ixdy  (3 .24) 


jAd 


exp 


jk_ 

2d 


((PAd)2+(qAd)2) 


3]  0(x,  y)  exp 


jk  2  ,  ..2  > 


2d 


tf+r) 


where  3(»)  indicates  the  forward  Fourier  transform.  The  substitutions  reveal  that/?  and  q 
are  the  frequency  domain  variables  of  x  and  y. 

£ 

Considering  the  £  dimension,  if  the  scan  is  limited  to  and  is  sampled# 

times,  then  the  sample  size  is 

which  in  the  p  dimension  is 


= 

A p  = 


L 


N 

A i 

Ad 


The  highest  spatial  frequency  that  can  be  resolved  is  then  [Strum,  Kirk,  1989,  p.  391] 

'N\_A4N_4^ 


cc  =  A/? 

For  L=  <fjmax ,  the  maximum  resolution  is 


2  Ad  2  Ad 


(3.25) 

(3.26) 

1] 

(3.27) 


,  1  2  Ad 

\=  — =  — . 
v  a  L 

where  the  subscript  Ip  is  used  to  indicate  the  resolution  is  in  terms  of  line  pairs. 
Therefore,  the  resolution  in  terms  of  lines  is 

h=K^mU 

2  2a  L 

For  example,  if  A=632.8nm,  <7=50mm,  and  L=0.5mm,  then 

,  Ad  (632.8nm)(150m/7z)  10r>0 

h  =  —  = - =  189.8wm . 

L  0.5 mm 

Thus,  the  highest  spatial  frequency  that  can  be  resolved  is 


1 


2  h  2(189.8  uni) 


-  =  2.634 


Ip 


mm 


(3.28) 

(3.29) 

(3.30) 

(3.31) 


3.2.2  Mosaicking 

From  the  preceding  development  on  the  maximum  resolution  obtainable  from  a 
digitized  portion  of  a  hologram,  it  is  clear  that  larger  scan  areas  are  needed  to  improve  the 
resolution  of  reconstructing  holograms.  One  approach  is  to  sample  overlapping  images 
from  the  hologram  and  compose  these  individual  images  into  a  larger  image.  This 
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process  is  commonly  referred  to  as  image  mosaicking.  Obviously,  accurate  placement  of 
the  individual  images  into  the  correct  location  is  needed  to  avoid  errors  in  the 
reconstruction. 

To  examine  the  impact  of  inaccurate  placement,  a  method  for  measuring  the 
effects  must  be  chosen.  There  are  a  number  of  methods  used  to  qualitatively  and 
quantitatively  measure  optical  systems  [Mouroulis,  Zhang,  1992].  One  of  the  more 
commonly  used  is  the  radius  of  encircled  energy  [Rivolta,  1986,  p.  2404].  This  is  the 
technique  chosen  for  evaluating  the  accuracy  needed  for  mosaicking  small  images 
together  to  form  a  larger  image. 

The  encircled  energy  of  radius  R  is  defined  as  the  fraction  of  the  total  energy  of 
the  point  spread  function  contained  within  the  radius  R  from  the  center  of  the  diffraction 
pattern  [Slepian,  1965,  p.  1111],  [Clark,  et.  al.,  1984].  A  commonly  used  figure  of  merit 
is  the  radius  that  encircles  84%  of  the  energy.  This  is  the  fraction  of  energy  contained 
within  the  Airy  disk  in  the  diffraction-limited  case  [Mouroulis,  Zhang,  1992].  However, 
as  will  be  seen  later,  selecting  a  radius  was  not  needed. 

Having  chosen  the  encircled  energy  as  the  method  to  analyze  mosaicking,  this 
method  must  now  be  applied  in  some  fashion  to  obtain  empirical  results.  Analysis  was 
performed  through  simulation.  The  program  code  written  to  perform  the  simulation  is 
provided  in  Appendix  A.  Comments  are  provided  within  the  code  to  help  in 
understanding  its  function.  A  description  of  the  approach  used  in  the  simulation  follows. 

The  approach  used  in  the  simulation  was  to  analyze  the  effects  of  mosaicking 
error  on  sine  waves.  A  perfect  mosaic  would  consist  of  several  segments  of  the  sine 
wave  matched  together  to  form  one  long  sine  wave.  However,  if  any  error  occurs  in  the 
mosaicking  process,  then  there  would  be  a  discontinuity  at  the  location  where  any  two 
mismatched  segments  were  joined  together.  Figure  3.6  illustrates  the  two  situations. 
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Mathematically,  each  discontinuity  can  be  modeled  as  a  change  in  phase.  Also, 
the  phase  change  at  each  successive  discontinuity  is  relatively  to  the  previous  phase 
change.  Suppose  the  second  segment  is  mismatched  with  respect  to  the  first  by  A0,  and 
that  the  third  segment  is  perfectly  matched  with  the  second  segment.  The  result  is  that 
the  third  segment  must  be  in  error  by  the  same  phase  amount  A0,  with  the  respect  to  the 
first  segment.  Therefore,  if  the  third  segment  is  mismatched  by  A (f>2  with  respect  to  the 
second,  then  its  mismatch  relative  to  the  first  segment  is  A0,  +  A02 . 


The  simulation  used  the  above  described  accumulation  of  phase  from  segment  to 
segment  to  represent  mosaicking  error  in  the  analysis.  Each  segment  was  given  a 
uniformly  random  phase  deviation  bound  by  +A (j)  with  respect  to  the  previous  segment. 
Thus,  the  actual  phase  deviation  of  each  segment  from  the  first  was  an  integer  multiple  of 
±A (f> .  For  instance,  the  second  segment  would  deviate  from  the  first  by  ±A(f>  and  the 
third  would  deviate  by  ±2A <j)  with  respect  to  the  first  segment,  and  so  on. 

Once  the  phase  for  each  segment  is  determined,  the  composite  signal  would  be 
generated.  Then,  a  fast  Fourier  transform  (FFT)  is  performed  on  the  composite  signal. 
The  result  is  then  squared  to  obtain  data  representing  energy.  The  form  of  this  result  is  a 
sine  function.  This  is  due  to  the  finite  size  of  the  signal.  Of  course,  if  the  signal  were 
infinite,  the  result  would  be  impulses  at  ±  the  frequency/of  the  sinewave.  Next,  the 
encircled  energy  is  computed  by  summing  the  values  at  increasingly  larger  distances 
from  the  ideal  frequency/.  Finally,  the  whole  process  is  repeated  many  times  and 
averaged  for  statistical  purposes. 


Data  from  the  simulations  are  graphed  in  Figs.  3.7  to  3.16.  There  are  6  curves  in 

n  %  3k 

each  graph.  The  maximum  phase  deviation  for  each  curve  is  0,  ±— ,  ±— ,  ± — ,  ±n, 

4  2  4 


± —  starting  with  the  top  curve.  The  dashed  curve  in  each  graph  is  the  ideal  encircled 
4 

energy  curve  for  one  segment.  If  mosaicking  is  to  provide  improvement,  then  the 
encircled  energy  curve  of  the  mosaicked  signal  should  lie  above  that  of  one  segment. 
From  these  graphs,  one  can  see  a  interesting  result.  As  the  allowed  phase  deviation  is 
increased,  the  curves  gradually  approach  the  ideal  one  segment  curve  and  coincide  with  it 
when  the  maximum  phase  deviation  is  ±n  .  This  is  true  in  every  case,  regardless  of  the 
segment  size,  frequency  of  the  signal,  or  number  of  segments.  As  the  maximum  phase 
deviation  is  increased  further,  the  curves  drop  below  the  ideal  curve  but  then  turn  back 
and  eventually  rise  above  the  ideal.  This  cycle  apparently  continues  for  even  larger  phase 
deviations. 


* 


U. 


The  implications  of  this  result  are  important.  The  worst  case  scenario  in  terms  of 
error  tolerance  is  when  the  signal  (or  image)  contains  frequencies  up  to  the  maximum 
allowed  of  half  the  sampling  frequency.  In  this  case,  one  pixel  would  represent  a  phase 
shift  of  k  .  Thus,  the  error  must  be  kept  within  one  pixel.  However,  the  n  phase  shift 
corresponds  to  an  encircled  energy  curve  equal  to  that  of  one  segment.  The  purpose  of 
mosaicking  is  to  improve  the  resolution  capability.  Therefore,  the  accuracy  must  be 
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Seg  L=1024,  #  of  seg=4,  f=0.066667,  iter=500,  maxdev=7.854, 0  pad  to  32768 
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Figure  3.9  Encircled  energy  plots  for  8  segments  of  1024  pixels  each  and  a  frequency  of 


500,  #  of  seg=2,  f=0.066667,  iter=500,  maxdev=7.854,  0  pad  to  32768 
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Figure  3,10  Encircled  energy  plots  for  2  segments  of  1500  pixels  each  and  a  frequency  of  1/15. 


Seg  L=1500,  #  of  seg=4,  f=0.066667,  iter=500,  maxdev=7.854, 0  pad  to  32768 
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Seg  L=1500,  #  of  seg=8,  f=0.066667,  iter=500,  maxdev=7.854, 0  pad  to  32768 
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Seg  L=1024,  #  of  seg=2,  f=0.49751,  iter=500,  maxdev=7.854,  0  pad  to  32768 
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Seg  L=1024,  #  of  seg=4,  f=0.49751,  iter=500,  maxdev=7.854,  0  pad  to  32768 
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Figure  3.15  Encircled  energy  plots  for  8  segments  of  1024  pixels  each  and  a  frequency  of  1/ 


66667,  iter=500,  maxdev=7.854, 0  pad  to  1048576 


n 


better  than  the  n  phase  shift.  If  —  is  selected  as  the  needed  phase  accuracy  in 

2 


placement,  then,  for  the  case  of  frequencies  at  the  Nyquist  limit,  this  means  that  the 
accuracy  must  be  better  than  one  pixel.  For  more  error  tolerance,  the  maximum 
frequency  in  the  image  must  be  even  less  than  the  Nyquist  limitation  of  half  the  sampling 
frequency. 


A  mosaicking  program  was  written  as  part  of  this  research  that  would  do  simple 
placement  of  images  utilizing  a  basic  correlation  technique  to  find  the  basic  match  among 
images.  Once  the  mosaicking  accuracy  requirement  became  clearer,  further  development 
of  this  program  was  abandoned.  The  code  for  this  program  is  provided  in  Appendix  B. 


3.2.3  Depth  of  Field 

Depth  of  field  is  the  amount  of  movement  of  the  object  plane  allowed  before  the 
stationary  image  plane  shows  deterioration  of  the  image.  Similarly,  depth  of  focus  is  the 
amount  of  movement  of  the  image  plane  allowed  before  the  stationary  appears  out  of 
focus  [Longhurst,  1967,  p.  306-308].  If  depth  of  focus  can  be  defined,  then  the  depth  of 
field  can  be  found  through  use  of  the  thin  lens  equation.  Figure  3.17  illustrates  the 
relationships  of  the  various  parameters. 


Figure  3.17  Depth  of  field  geometry. 


The  depth  of  focus  in  an  optical  system  is  closely  related  to  its  resolution  ability. 
Just  as  the  Rayleigh  criterion  provides  a  measure  of  resolution,  there  is  a  similar  measure 
along  the  axis  of  the  optical  system.  Essentially,  the  light  distribution  at  the  focal  point 
must  be  considered  in  three  dimensions  instead  of  two  [Bom,  Wolf,  1980,  p.  435].  The 
intensity  along  the  axis  for  a  circular  aperture  is  of  the  form  [Bom,  Wolf,  1980,  p.  441] 

sinu/4 
m/4 


(3.32) 


where  [Bom,  Wolf,  1980,  p.  437] 


(3.33) 


In  Eq.  (3.33),  z'  is  the  optical  axis,  r  is  the  radius  of  the  aperture,  and  R  is  the  radius  of 
the  wavefront  at  the  aperture  [Born,  Wolf,  1980,  p.  435-436].  An  intensity  loss  of  20% 
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of  that  at  the  center  is  considered  acceptable  [Bom,  Wolf,  1980,  p.  441].  This 
corresponds  to 

ms  3.2.  (3.34) 

Substituting  Eq.  (3.34)  into  Eq.  (3.33)  gives  the  depth  of  focus  as  [Bom,  Wolf,  1980,  p. 
441] 

1  (  R\2 

As'  =  ± - A.  (3.35) 

2\rJ 

This  is  usually  combined  with  the  Rayleigh  criterion  of  resolution,  which  can  be  stated  as 
[Hecht,  1998,  p.  464] 

1.22/A 


K=- 


2  r 


to  give  the  depth  of  focus  as 


*  ,  2.69 V 
As  = — -1 - 


(3.36) 

(3.37) 


where  Hr  is  the  Rayleigh  criterion  resolution  of  the  system.  Note  that  the  radius  of 
curvature  R  in  Eq.  (3.33)  is  estimated  to  be  equal  the  focal  length /in  Eq.  (3.36).  The 
depth  of  focus  is  commonly  approximated  as  [Vikram,  1992,  p.  62] 


A 


(3.38) 


The  above  development,  as  mentioned,  is  based  on  the  Rayleigh  criterion 
resolution  limit  of  a  system.  However,  from  the  development  earlier  on  the  maximum 
theoretical  resolution  due  to  a  limited  scan  area,  the  actual  resolution  is  worse  than  the 
Rayleigh  criterion.  Therefore,  it  is  more  appropriate  to  use  Eq.  (3.35)  directly  for  our  non 
diffraction-limited  case.  If  R  in  Eq.  (3.35)  is  assumed  to  be  roughly  equal  to  the  distance 
d,  and  assuming  2r  is  equal  to  the  length  L  of  the  limited  scan  area,  then  the  resulting 
depth  of  focus  is 

Ay/  =  ±2A^_  (3.39) 


All  of  the  developments  above  make  the  assumption  that  the  aperture  is  circular. 
However,  since  a  video  camera  is  being  used  to  digitize  the  hologram,  the  aperture  is 
actually  rectangular  because  the  electronic  imager  is  inherently  rectangular.  Thus,  a 
development  similar  to  Bom  and  Wolf  cited  earlier  except  based  on  a  rectangular 
aperture  should  be  applied.  Indeed,  Bom,  Wolf,  1980,  p.  440  cited  Matthews  and  Cullen 
for  their  study  of  a  square  microwave  lens.  Their  work  for  a  uniformly  illuminated 
square  aperture  can  be  summarized  as  follows  [Matthews,  Cullen,  1956,  p.  449-453]. 
The  scalar  wave  function  was  determined  to  be  [Matthews,  Cullen,  1956,  p.  452-453] 


Ad  q 


where  C(x)  and  S(x)  are  Fresnel  integrals  defined  by 


C(x)  +  jS(x)  =  exp 


(3.40) 


(3.41) 


35 


and 


2 


q  d 2 


(3.42) 


The  magnitude  and  phase  of  Eq.  (3.40)  applies  for  q  positive,  however  the  magnitude  is 
correct  for  q  positive  or  negative. 

The  depth  of  focus  is  derived,  as  in  Bom  and  Wolf,  by  considering  the  drop  in 
amplitude  in  the  z'  direction.  The  first  step  is  to  determine  the  intensity  of  the  wave  at 
z'  =  0  from  Eq.  (3.40).  This  was  determined  by  substituting  Eqs.  (3.41)  and  (3.42)  into 

Eq.  (3.40)  and  calculating  |0(O)|2  which  gave  (with  the  help  of  Mathcad  by  Mathsoft) 


tp(-yfe') 


2  z%7 


(3.43) 


16r?; 

Using  the  same  criterion  for  depth  of  focus,  that  is  when  the  intensity  has  dropped  to  0.8 
of  the  intensity  at  the  focal  plane,  then  the  edges  of  the  depth  of  focus  will  be  located 
where  z  satisfies 


I 


(3.44) 


(3.45) 


Unfortunately,  a  nice  closed  form  of  Eq.  (3.44)  has  not  been  found.  Therefore,  it 
is  necessary  solve  Eq.  (3.44)  through  a  numerical  technique  which,  for  example,  finds  the 
root  of 

=  (3-45) 

As  an  example,  let  X=632.8nm,  Tj0^466.3fJm,  and  d=73.25mm.  (Note  that  2rj0  is  the 
dimension  of  one  side  of  the  aperture.)  Substituting  these  into  Eq.  (3.45)  and  solving  for 
z'  (with  the  help  of  Mathcad)  gives  z'  =  5547mm .  Thus,  the  depth  of  focus  is 

As'  =  ±5547 mm .  (3.46) 

In  comparison,  if  the  aperture  was  circular  with  a  diameter  of  932.6pm  (the  same  size  as 
one  side  of  the  square  aperture),  then  the  depth  of  focus  given  by  Eq.  (3.39)  is 

As'  =  ±7.808mm.  (3.47) 

To  determine  the  relation  between  depth  of  focus  and  depth  of  field,  we  start  with 
the  thin  lens  equation,  which  can  be  written  as 

(3.48) 

s  s  f 

where  s  is  the  object  distance  from  the  lens,  s'  is  the  image  distance,  and/is  the  focal 
length  of  the  lens.  Equation  (3.48)  can  be  rewritten  as 

<3-49> 

s  f  s 
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The  derivative  of  Eq.  (3.49)  then  gives  the  relation  between  depth  of  field  and  depth  of 
focus  as 


ds  =  -—ds'. 
s 


(3.50) 


For  objects  far  away  from  the  lens,  then  s'  is  approximately  equal  to  the  focal  length  of 
the  lens,  and  so  Eq.  (3.50)  can  be  written  as  approximately 

ds  =  — —rds' .  (3.51) 

f2 


To  determine  the  depth  of  field  in  holograms,  it  is  helpful  to  depict  the  reconstruction  of 
the  real  image  from  a  hologram  as  shown  in  Fig.  3.18  [Bom,  Wolf,  1980,  p.  454],  [Kreis, 
1996,  p.  157].  In  the  case  of  digital  reconstruction,  the  lens  can  be  numerically  simulated 
and  positioned  exactly  at  the  hologram  plane  [Kreis,  1996,  p.  157].  If  no  magnification  is 
desired,  then  s=  s'  and  so  the  magnitudes  of  depth  of  field  and  depth  of  focus  are  equal. 
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4.0  EXPERIMENTAL  DESIGN 


4.1  Objective  of  Experiments 

There  are  two  main  objectives  to  be  accomplished  through  experimentation.  The 
first  is  to  verify  the  resolution  achievable  for  a  limited  digitized  area  of  a  hologram,  as 
discussed  in  Chapter  3.  The  approach  will  be  to  make  a  hologram  of  a  resolution  target. 
Such  targets  have  lines  or  bars  of  known  spatial  resolution.  The  second  goal  is  to 
demonstrate  the  depth  of  field  expected  in  digital  reconstructions.  The  approach  will  be 
to  make  a  hologram  of  a  custom  designed  “ruler”  that  is  oriented  at  an  angle  with  respect 
to  the  hologram  plane.  In  this  way,  certain  parts  of  the  ruler  depending  on  the 
reconstruction  distance. 

4.2  Components 

Before  the  experiments  can  be  designed,  the  hardware  available  for  use  must  be 
described  since  their  specifications  will  be  a  factor  in  the  actual  designs.  The  object  used 
to  measure  resolution  is  a  standard  United  States  Air  Force  1951  target.  A  close-up  view 
of  the  target  is  provided  in  Fig.  4. 1 .  The  chart  has  sets  of  three  lines,  both  horizontal  and 
vertical,  that  range  in  spatial  frequency  from  1.00  to  228.0  line  pairs  per  millimeter 
(lp/mm).  The  camera  used  for  digitizing  the  holograms  is  a  Kodak  Megaplus  1.4  which 
uses  a  CCD  imager  having  1316  x  1033  pixels  where  each  pixel  is  6.8  um  square  and  the 
fill  ratio  is  100  percent.  A  Matrox  frame  grabber  was  used  for  obtaining  and  storing 
digitized  images  from  the  camera.  A  standard  Pentium-based  personal  computer  was 
used  communicate  with  the  frame  grabber  and  store  the  digitized  images  onto  disk.  A 
Nikon  microscope  with  a  10X  objective  having  a  numerical  aperture  of  0.25  was  used. 
The  magnification  M  was  determined  experimentally  to  be  7.4664.  How  this  was 
determined  will  be  discussed  in  the  next  chapter.  A  Hughes  model  3224H-PC  helium- 
neon  laser  was  used  to  make  the  holograms  and  so  the  wavelength,  X ,  is  632.8nm. 


Figure  4.1  United  States  Air  Force  1951  resolution  target. 
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4.3  Resolution  Test 


The  objective  of  this  experiment  is  to  demonstrate  that  the  resolution  limitation 
due  to  a  limited  digitized  area  of  a  hologram  is  as  described  by  Eq.  (3.29).  The  approach 
to  verify  this  equation  through  experimentation  is  illustrated  in  Fig.  4.2.  There  are  several 
parameters  of  the  experimental  setup  that  must  be  taken  into  consideration  to  successfully 
demonstrate  the  resolution  achievable.  These  parameters  are  the  distance  of  the 
resolution  chart  from  the  film  plate,  the  frequency  of  the  resolution  bars  on  the  chart,  the 
orientation  and  wavelength  of  the  object  and  plane  waves,  the  numerical  aperture  and 
magnification  of  the  microscope  objective,  and  the  wavelength  of  the  microscope 
illumination.  Explanation  of  each  of  these  parameters  now  follows. 


chart 


Figure  4.2  Resolution  test  setup. 


The  first  parameter  considered  is  the  frequency  of  the  resolution  bars  on  the  chart. 
The  sets  of  lines  that  can  be  resolved  in  this  experiment  are  determined  by  Eq.  (3.29). 
Since  a  microscope  is  used  to  digitize  the  hologram,  the  magnification  must  be  included 
in  calculating  L  in  Eq.  (3.29).  Defining  C  as  the  actual  size  of  the  imager,  then 


(4.1) 


The  number  of  pixels  in  the  imager  multiplied  by  their  size  gives  C.  Substituting  Eq. 

(4.1)  into  Eq.  (3.29)  gives 


C 

~  2  MAd  ’ 


(4.2) 
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where  the  subscript  o  indicates  the  frequency  is  for  the  object  (as  opposed  to  the 
hologram  which  will  be  addressed  later). 

Substituting  all  of  the  parameters  given  earlier  in  Eq.  (4.2)  and  taking  into 
consideration  the  practicality  of  the  distance  the  target  could  be  placed  from  the  film 
plate,  it  was  decided  to  place  the  target  at  distances  of  164mm,  218mm,  and  350mm. 
Actually,  the  hologram  of  the  target  at  distances  of  164mm  and  350mm  was  done  as  a 
double  exposure  for  the  purpose  of  the  depth  of  field  experiment.  The  expected 
resolution  of  the  target  achievable  for  the  distances  of  164mm,  218mm,  and  350mm  were 
4.493  lp/mm,  3.380  lp/mm,  and  2.105  lp/mm,  respectively.  On  the  resolution  target, 
these  spatial  frequencies  correspond  to  Group  2  Element  2,  Group  1  Element  5,  and 
Group  1  Element  1,  respectively. 


For  this  same  experiment,  there  is  another  spatial  frequency  resolution  problem 
that  must  be  addressed.  That  is  the  spatial  frequency  of  the  fringes  created  on  the 
hologram  due  to  the  angles  of  the  object  and  reference  waves.  The  spatial  frequency 
resulting  from  the  interference  of  the  object  and  reference  waves  was  derived  in  Chapter  2 
and  the  result  [Eq.  (2.34)]  was 


«*  = 


_  sin(flfl)  +  sin(fl„) 
X 


(4.3) 


where  the  subscript  h  is  used  to  signify  the  hologram.  From  Fig.  4. 1 ,  the  geometry  of  the 
setup  can  be  used  to  derive  the  angles  6a  and  9U  as 


e 


a 


=  tan 


f-1 

\d) 


(4.4) 


and 


6„  =  tan 


-i 


x— o') 

~d~J 


Substituting  Eqs.  (4.4)  and  (4.5)  into  Eq.  (4.3)  gives 


sin 

«*=  — 


tan 


(3 


f 


+  sin 


tan 


-i 


V 


(t)) 


(4.5) 


(4.6) 


The  importance  of  Eq.  (4.6)  is  that  a  microscope  objective  has  a  limit  to  its 
resolution  capability  usually  expressed  as  the  numerical  aperture  (NA)  as  discussed  in 
Chapter  3.  Therefore,  the  geometry  of  the  experimental  setup  must  be  adjusted  so  that  the 
microscope  can  resolve  the  interference  fringes  recorded  on  the  hologram.  The  NA 
required  to  resolve  the  spatial  frequencies  described  by  Eq.  (4.6)  can  be  expressed  by 
combining  Eqs.  (3.7)  and  (3.29)  which  gives 

NA  =  0.6l-2Xm-ah,  (4.7) 

where  the  subscript  m  is  used  to  indicate  that  Xm  is  the  wavelength  of  the  light  used  by 
the  microscope  to  illuminate  the  hologram.  The  microscope  contains  a  green  interference 
filter  with  a  specified  center  frequency  of  546nm.  The  filter  bandwidth  around  this  center 
frequency  is  not  known.  The  longest  wavelength  in  the  bandwidth  would  correspond  to 
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the  most  stringent  NA  requirement.  The  longest  wavelength  is  estimated  to  be  600nm, 
and  so  that  was  the  wavelength  used  for  Xm  in  Eq.  (4.7). 

Table  4.1  provides  a  number  of  NAs  for  various  object  and  hologram  locations  to 
show  the  variations  that  exist.  The  first  four  parameters  in  Table  4.1,  namely  d,  v,  x,  and 
o,  are  shown  in  Fig  4.2.  The  parameter  d  is  the  distance  the  resolution  target  is  from  the 
hologram  film  plate.  The  parameter  v  is  the  distance  the  edge  of  the  resolution  target  is 
from  the  edge  of  the  hologram  film  plate.  The  parameters  x  and  o  are  specific  locations 
on  the  hologram  plate  and  resolution  target,  respectively.  As  can  be  seen,  the  NA  can 
vary  significantly,  depending  on  what  part  of  the  hologram  is  being  digitized.  Because  of 
these  variations,  it  is  important  to  consider  the  geometry  of  the  setup  when  deciding  what 
part  of  the  hologram  to  digitize. 


Table  4.1  Effect  of  various  parameters  on  the  numerical  aperture. 


d  (mm) 

o(mm) ' 

X  (nm)  °h  i 

■  Ip/inm)  i 

V.  .  w 

W(nm) 

NA  . 

164.00 

34.00 

0.00 

0.00 

632.80 

320.80 

600.00 

0.23 

164.00 

34.00 

0.00 

15.00 

632.80 

176.86 

600.00 

0.13 

164.00 

34.00 

0.00 

34.00 

632.80 

0.00 

600.00 

0.00 

164.00 

34.00 

15.00 

0.00 

632.80 

464.73 

600.00 

0.34 

164.00 

34.00 

15.00 

15.00 

632.80 

320.80 

600.00 

0.23 

164.00 

34.00 

15.00 

34.00 

632.80 

138.93 

600.00 

0.10 

164.00 

34.00 

34.00 

0.00 

632.80 

641 .59 

600.00 

0.47 

164.00 

34.00 

34.00 

15.00 

632.80 

502.66 

600.00 

0.37 

164.00 

34.00 

34.00 

34.00 

632.80 

320.80 

600.00 

0.23 

Figure  4.3  is  a  picture  of  the  actual  setup  used  in  the  experiment.  In  order  to 
facilitate  the  positioning  of  the  target,  hologram  film  plate,  object  wave,  and  plane  wave, 
the  edges  of  the  target  and  hologram  film  plate  were  used  as  guides.  This  can  be  seen  by 
referring  back  to  Fig.  4.2.  In  Fig.  4.2,  the  long  vector  of  the  reference  plane  wave 
illustrates  how  these  edges  were  used.  The  reference  plane  wave  was  positioned  so  that  it 
just  barely  missed  the  edge  of  the  resolution  target  and  landed  on  the  edge  of  the 
hologram  as  shown.  This  made  positioning  the  plane  wave  at  the  desired  angle  easier. 

The  entire  resolution  target  was  not  placed  directly  in  front  of  the  film  plate.  The 
target  and  film  holder  were  adjusted  so  that  the  portion  of  the  target  with  the  lines  just 
appeared  at  the  edge  of  the  film  plate.  This  allowed  using  a  smaller  plane  wave  angle  to 
reduce  the  spatial  frequencies  on  the  hologram. 

4.4  Depth  of  Field  Test 

In  order  to  test  the  depth  of  field,  two  different  approaches  were  taken.  First,  a 
double  exposure  hologram  was  made  with  the  resolution  target  placed  at  two  different 


Figure  4.3  Resolution  test  setup. 


distances  from  the  hologram  film  plate.  The  second  approach  was  to  make  a  hologram  of 
a  ruler  at  a  slanted  angle.  The  depth  of  field  can  be  calculated  by  solving  for  the  z' 

which  satisfies  Eq.  (3.44)  which  when  combined  with  Eq.  (4.1)  and  noting  that  r\0  -  — 


becomes 


The  depth  of  field  is  then 


:  0.8  C 

A2  mV4 

As  =  ±z' . 


=  0. 


(4.8) 

(4.9) 


For  the  first  approach,  the  setup  is  the  same  as  in  Figs.  4.2  and  4.3  except  that  the 
resolution  target  was  placed  at  two  different  locations.  Since  depth  of  field  is  more  of  a 
subjective  parameter,  it  is  difficult  to  measure.  So,  the  two  distances  chosen  were  such 
that  the  degradation  would  definitely  be  noticeable.  Since,  a  hologram  with  distance  of 
164mm  from  the  hologram  film  plate  had  already  been  chosen  for  the  resolution  test,  it 
was  decided  to  use  this  as  one  of  the  distances.  At  this  distance,  the  total  depth  of  field 
from  Eqs.  (4.8)  and  (4.9)  is 

As  =  ±27 ,8mm.  (4.10) 

This  indicated  that  the  second  distance  would  have  to  be  at  least 

1 64mm + 27.8mm  =  1 9 1.8mm .  (4.11) 
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However,  since  depth  of  field  increases  as  the  distance  increases,  then  the  distance  would 
have  to  be  even  further.  Another  consideration  is  that  the  resolution  target  would  need  to 
be  resolvable  at  the  second  distance.  Based  on  this,  a  distance  of  350mm  was  selected. 

At  this  distance,  Group  1  Element  1  should  be  resolvable.  The  corresponding  depth  of 
field  would  be 

Ay  =  ±1 26.6 mm .  (4.12) 

Figure  4.4  shows  the  relation  between  the  two  depths  of  field  and  distances. 


|«-As=55.6  lmm-i ►] 
As=253.3mm - ►! 


H —  d=  164mm  — ► 
« -  d=350mm - ► 


Figure  4.4  Illustration  of  two  different  depths  of  field. 

The  second  approach  to  testing  the  depth  of  field,  as  mentioned  earlier,  was  to  use 
a  ruler  at  a  slanted  angle.  Figure  4.5  illustrates  this  experiment.  One  limitation  that  had 
to  be  taken  into  account  was  the  resolution  issue.  The  size  of  the  ruler  markings  had  to  be 
such  that  they  could  be  resolved  at  the  furthest  distance  the  ruler  would  be  placed  from 


Figure  4.5  Depth  of  field  test  setup. 
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the  hologram  film  plate.  This  ruled  out  use  of  traditional  rulers.  So,  a  ruler  was  designed 
and  printed  on  a  transparency  using  a  laser  printer.  The  ruler  used  is  shown  in  Fig.  4.6. 
Note  that  three  of  the  lines  are  only  half  the  length  of  the  others.  This  was  done  to  help  in 
identifying  what  part  of  the  ruler  has  been  reconstructed. 


* 


4 


4 


Laser 

— ► 


Figure  4.6  Line  pattern  used  as  a  ruler  in  one  depth  of  field  test. 


For  this  experiment,  there  were  a  number  of  parameters  to  consider  in  the  design 
process.  Besides  the  NA,  it  was  also  necessary  to  consider  the  angle  of  the  ruler  and  how 
that  would  affect  the  apparent  size  of  the  individual  ruler  lines.  The  angle  was  somewhat 
controlled  by  the  width  of  the  object  plane  wave.  The  NA  affected  the  width  and  the 
length  of  the  ruler.  After  consideration  of  all  of  these  factors,  it  was  decided  that  the 
distance  would  be  set  to  185mm,  the  width  would  be  set  to  24mm,  and  the  length  of  the 
ruler  would  be  set  to  1 14.3mm. 


Determining  the  NA  for  this  setup  is  not  as  straightforward  as  in  the  resolution 
target  setup.  This  is  because  the  object  (ruler)  varies  in  distance  from  the  hologram  film 
plate.  From  Fig.  4.5,  the  angle  of  the  reference  wave  is  given  by 


0a  =  a  tan| 


(4.13) 


The  angle  of  the  object  wave  is  given  by 


6„  =  atan 


x-o 


(v-o) 


d- 


tan(0r) 


(4.14) 


Table  4.2  shows  how  the  NA  varies  for  this  geometry.  Figure  4.7  is  a  picture  of  the  setup 
used  for  this  depth  of  field  test.  Similar  to  the  resolution  test,  the  reference  plane  wave 
was  adjusted  so  that  it  just  barely  missed  the  edge  of  the  ruler  farthest  from  the  film  plate. 
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Table  4.2  Numerical  aperture  values  for  depth  of  field  test  setup. 


;  d  (ram) . 

v  (mm) 

x  (mm)  <*  (inm) 

••X  (tmi) 

a*,  (ip/mm) 

Ha  (nn»1| 

HP! 

185.00 

24.00 

0.00 

0.00 

632.80 

203.31 

600.00 

0.15 

185.00 

24.00 

0.00 

15.00 

632.80 

38.55 

600.00 

0.03 

185.00 

24.00 

0.00 

24.00 

632.80 

0.00 

600.00 

0.00 

185.00 

24.00 

15.00 

0.00 

632.80 

520.34 

600.00 

0.38 

185.00 

24.00 

15.00 

15.00 

632.80 

203.31 

600.00 

0.15 

185.00 

24.00 

15.00 

24.00 

632.80 

126.52 

600.00 

0.09 

185.00 

24.00 

24.00 

0.00 

632.80 

695.34 

600.00 

0.51 

185.00 

24.00 

24.00 

15.00 

632.80 

302.50 

600.00 

0.22 

185.00 

24.00 

24.00 

24.00 

632.80 

203.31 

600.00 

0.15 

Film  plate 
holder 


“Ruler”  and 
scatter  plate 


Beam  splitter 


Aperture 


5.0  DIGITAL  HOLOGRAPHY  RECONSTRUCTION  PROGRAM 


In  this  chapter,  the  mathematics  and  software  code  used  in  the  digital 
reconstruction  process  will  be  explained.  The  software  used  in  the  development  of  this 
program  is  called  IDL  by  Research  Systems,  Inc.  In  this  section  only  the  essential  math  in 
the  process  will  be  addressed.  The  full  program  including  all  of  the  data  I/O,  GUIs,  etc  is 
included  in  Appendix  C.  In  Chapter  2,  the  equation  for  calculating  the  wavefront  at  the 
hologram  plane  was  given  as 

M(4  +  ”  >J 


U(^t?)  =  TTTexP 


jAd 
JJ  0(x,  y)  exp 


exp 

L  / 


j2n 

Ad 


(x^  +  yrj)  ]dxdy 


(5.1) 


In  the  reconstruction  process,  Eq.  (5.1)  can  also  be  used  except  that  now  the  equation 
becomes 


Jkd 


0(x',/)  =  — exp| 


•JJU(£,i7)exp 


jja 

2d 


eXP 


-jjj-ix'Z+y'ri)  d&ri. 


(5.2) 


where  U(<i;,77)  is  the  wavefront  at  the  hologram  plane  and  0(x',y')  is  the  wavefront  at 
the  reconstruction  plane.  In  order  to  obtain  U(<^,tj)  for  use  in  the  calculations,  the 
digitized  portion  of  the  hologram  must  be  multiplied  by  the  reference  plane  wave.  The 
result  is  U(£,77) .  Next  in  the  equation  is 

jk 


exp 


Id 


(f  +*f) 


(5.3) 


Equation  (5.3)  is  sometimes  referred  to  as  a  Gaussian  phase  term.  Multiplying  Eq.  (5.3) 
with  U(£,77)  and  then  taking  the  forward  Fourier  transform  gives  the  double  integral 
portion  of  Eq.  (5.1).  Multiplying  this  by  the  terms  in  front  of  the  double  integrals  in 
Eq.(5.1)  gives  the  final  desired  result  0(x',y') . 


Mathematically  creating  a  reference  plane  wave  begins  with  the  basic  equation  for 
describing  electromagnetic  (EM)  wave  propagation.  Defining  the  reference  plane  wave 
as  P(^,?7) ,  the  general  equation  is 

P(£,?7) =  exp  j(k  -r-wt  +  fj)) .  (5.4) 


The  last  two  factors  in  the  exponent  can  be  ignored.  The  dot  product  is  given  by 
expressing  k  in  its  £  and  77  components.  Using  the  spherical  coordinate  system,  the 
plane  wave  can  be  described  by 


P(£,  77)  =  exp  j[k  cos  (6)  sin  (<j>)  +  k  sin(0)  sin(0)] . 

(5.5) 

Defining 

k^  =  k  cos(0)  sin  (<p) 

(5.6) 

and 
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Eq.  (5.5)  becomes 


kn=k  sin(0)  sin  (0) , 

P(£,77)  =  exp  7(*{+V- 


(5.7) 


(5.8) 


The  software  code  that  creates  the  plane  wave  is  the  following  (note  that 
throughout  the  software  code  the  letters  E,  and  77  were  replaced  with  x  and  y  because  of  the 
difficulty  in  using  Greek  letters  in  the  code): 


k=2.*  !dpi/info.lambda 

kx=k*cos(info. theta*  !dpi/180d)*sin(info.phi*  !dpi/l  80d) 
ky=k*sin(info.theta*!dpi/180d)*sin(info.phi*!dpi/180d) 
vectx=(dindgen(info.mc)-round(info.mc/2))*info.dx/info.mag 
arx=transpose(make_array(info.nc,  1 , /double,  value=l  .0)) 
vecty=(transpose(dindgen(info.nc)-round(info.nc/2)))*info.dy/info.mag 
ary=make_array(info.mc,  1 , /double,  value=  1 .0) 
planewave=kx  *  (vectx#arx)+ky  *  (ary#vecty) 
planewave=exp(dcomplex(0,temporary(planewave))) 


In  the  above  code,  the  variable  info  is  a  named  structure  that  is  basically  a  container  of 
variables.  So,  info.lambda  refers  to  the  variable  lambda  contained  in  the  structure  info. 
The  variables  in  info  called  in  the  above  code  are  lambda  which  is  the  laser  wavelength, 
theta  and  phi  which  define  the  orientation  of  the  plane  wave  in  spherical  coordinates,  me 
and  nc  which  are  the  dimensions  of  the  digitized  hologram,  dx  and  dy  which  are  the  actual 
size  of  the  pixels  in  the  camera,  and  finally  mag  which  is  the  actual  magnification  of  the 
microscope.  In  addition,  !dpi  calls  the  double  precision  value  of  K .  Whenever  d  is  used 
next  to  a  number  as  in  180d,  this  defines  the  number  as  double  precision  floating  point 
number.  The  result  of  the  code  is  that  it  calculates  the  plane  wave  as  illustrated  in  the 
equation  where  the  digitized  hologram  size  is  assumed  to  be  1024  X  1024  pixels: 


planewave  =  exp  j 
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(5.9) 


Multiplying  the  result  of  Eq.  (5.9)  with  the  digitized  portion  of  the  hologram  creates 
U(£,77)  in  Eq.  (5.2). 


The  next  term  in  Eq.  (5.2)  after  U(^,q)  is  the  Gaussian  phase  term  given  in  Eq. 
(5.3).  If  the  original  plane  wave  and  not  the  conjugate  plane  wave  is  used  to  recreate  the 
object  wavefront  at  the  hologram  plane,  then  a  Fourier  lens  must  be  introduced  to  focus 
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the  real  image  at  the  observation  plane.  The  phase  transformation  of  a  lens  with  a  focal 
length  of/is  [Goodman,  1996,  p.  99] 


=  exp 


(5.10) 


Taking  the  approach  of  multiplying  the  wavefront  at  the  hologram  plane  by  the 
transformation  of  Eq.  (5.10)  results  in  essentially  placing  a  lens  of  focal  length /  at  a 
distance  d  from  both  the  object  and  image.  Thus,  the  focal  length  of  the  lens  is 

/=  f.  (5.11) 

Mathematically,  the  lens  transformation  Eq.  (5.10)  is  accounted  for  by  placing  it  into  Eq. 
(5.2)  resulting  in  [Goodman,  1996,  p.  102-104] 
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Substituting  Eqs.  (5.10)  and  (5.1 1)  into  Eq.  (5.12)  gives 
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(5.13) 


Note  that  the  difference  between  Eqs.  (5.1)  and  (5.13)  is  the  Gaussian  phase  term  of  Eq. 
(5.1)  is  replaced  by  its  conjugate  to  give  Eq.  (5.13). 


The  Gaussian  phase  term  of  Eq.  (5.13)  is  created  by  the  following  software  code: 


gauss=-((vectx#arx)A2+(ary#vecty)A2)*k/2./info.d 

gauss=exp(complex(0,temporary(gauss))). 


The  variable  info.d  refers  to  the  distance  of  the  object  from  the  hologram  (and  also  the 
distance  from  the  hologram  to  the  real  image).  The  Gaussian  phase  term  as  calculated  by 
the  software  code  can  be  illustrated  as  follows  where  again  the  digitized  hologram  size  is 
assumed  to  be  1024  X  1024  pixels: 


gauss  =  exp 
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The  final  steps  in  reconstructing  the  image  from  the  hologram  are  fairly 
straightforward.  The  software  code  lines  are  as  follows: 

recon=double(holo)*gauss*planewave 
recon=fft(temporary(recon),- 1  ,double=  l , /overwrite) 
recon=abs(temporary(recon))A2 
recon=shift(temporary(recon),info.mc/2,info.nc/2). 

The  first  step  is  to  multiply  the  digitized  portion  of  the  hologram  (called  hold),  the  plane 
wave  calculated  in  Eq.  (5.9),  and  the  Gaussian  phase  term  in  Eq.  (5.14)  together.  The 
result  is  called  recon.  Then  a  forward  Fourier  transform  of  recon  is  calculated  and  placed 
back  into  recon  (to  save  memory).  The  absolute  value  squared  of  recon  is  calculated 
since  the  desired  result  is  the  intensity  because  that  is  what  we  see  with  our  eyes.  Finally, 
the  image  is  shifted  because  the  discrete  Fourier  transform  gives  the  positive  frequencies 
first  followed  by  the  mirror  image  or  negative  frequencies. 

The  Graphical  User  Interface  (GUI)  is  shown  in  Fig.  5.1.  The  general  process  is 
to  load  the  interference  image  to  be  reconstructed  and  then  input  the  various  parameters 
such  as  the  distance  of  the  object  from  the  hologram  film  plate,  magnification  of  the 


Figure  5.1  Graphical  User  Interface  (GUI)  of  the  reconstruction  program. 


microscope,  size  of  the  pixels,  etc..  The  buttons  for  partial  images  allow  cutting  the  size 
of  an  image  so  that,  for  instance,  the  size  can  be  adjusted  to  powers  of  2.  The  Linear/Log 
Display  button  is  used  to  switch  between  the  raw  reconstructed  image  being  displayed  or 
to  take  the  log  of  the  raw  reconstructed  image  and  then  displaying  it.  The  latter  helps  to 
bring  out  detail  in  reconstructions  that  have  a  large  dynamic  range.  The  Mask  Image 
Portion  provides  an  additional  method  of  bringing  out  detail  in  reconstructed  images.  It 
allows  brighter  parts  of  the  image  to  be  suppressed  by  essentially  drawing  a  black  box  on 
top  of  them. 


* 


* 
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6.0  EXPERIMENTAL  RESULTS 

6.1  Microscope  Magnification 

As  discussed  in  Chapter  4,  the  magnification  of  the  microscope  objective  must  be 
included  in  the  calculations  because  it  directly  affects  the  size  of  the  hologram  digitized 
by  the  camera.  Therefore,  the  actual  amount  of  magnification  must  be  determined.  The 
method  of  determining  the  magnification  was  to  digitize  the  resolution  target  itself  with 
the  microscope  and  camera.  Since  the  sizes  of  the  lines  on  the  resolution  target  are 
known,  then  the  magnification  can  be  determined.  In  order  to  reduce  the  error  as  much  as 
possible,  the  largest  resolution  target  line  that  would  fit  into  the  field  of  view  was  chosen. 
This  was  determined  to  be  Group  0  Element  1,  this  largest  on  the  chart.  Figure  6.1  shows 
the  digitized  image.  The  brighter  pixels  are  the  line,  which  fills  about  half  of  the  vertical 
field  of  view. 


Figure  6.1  Digitized  image  of  a  Group  0  Element  1  line  from  the  resolution  target. 

* 

From  the  specifications  of  the  resolution  target,  it  is  known  that  Group  0  Element 
1  lines  are  1 .00  line  pairs/mm.  Therefore,  one  line  is  0.5mm  in  width.  If  no 
magnification  were  used,  then  one  line  would  correspond  to  4 

fOSmnY lpixd\i5ipixeh  (dl) 

V  line  yv  6.8 um )  line 

In  other  words,  the  resolution  target  line  would  subtend  73.53  pixels  on  the  camera’s 
imager.  Now,  if  the  number  of  pixels  subtended  by  the  same  line  under  the  microscope 
can  be  determined,  .then  the  magnification  can  be  calculated.  Using  MATLAB  to  read 
and  zoom  in  on  Fig.  6.1,  the  number  of  pixels  can  be  determined.  Figures  6.2  and  6.3 
show  the  upper  and  lower  edges,  respectively.  At  this  time,  the  reader  is  referred  to 
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Appendix  D  which  discusses  relative  error  propagation.  It  is  necessary  in  understanding 
how  the  errors  discussed  in  the  remainder  of  this  chapter  are  estimated  and  calculated. 

Figures  6.2  and  6.3  illustrate  the  difficulty  in  determining  the  exact  magnification. 
Deciding  the  location  of  the  edges  precisely  is  not  possible.  This  is  where  the  error 
analysis  and  propagation  discussed  earlier  must  be  applied.  From  Fig.  6.2,  it  appears  that 
the  upper  edge  is  located  at  pixel  number  368.  However,  this  may  not  be  the  true  edge. 

A  reasonable  estimate  of  the  error  is  ±  2  pixels.  From  Fig.  6.3,  it  appears  that  the  lower 
edge  is  located  at  pixel  number  917.  Again,  the  estimate  of  error  is  ±  2  pixels.  The 
resulting  width  and  corresponding  relative  error  of  the  line  is  then 

549 pixels+Apixels  =  550pixels±0.12S6% .  (6.2) 

(Keep  in  mind  that  the  relative  error  is  actually  approximated  since  the  measured  as 
opposed  to  the  true  value  is  used  as  discussed  earlier  in  this  chapter.)  The  actual 
magnification  of  the  microscope  objective  can  now  be  calculated  to  be 

M  =  549P**e/5± 0.7286%  _  ?  4664  ±  o  7286%  .  (6.3) 

15 .55  pixels 
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Figure  6.2  Magnified  upper  edge  of  resolution  line. 

6.2  Unfeasibility  of  Testing  Mosaicking 

The  significant  result  from  section  3.2.2  on  mosaicking  is  that  accuracy  of 
placement  must  be  better  than  one  pixel  to  achieve  an  improvement  in  resolution.  One 
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Figure  6.3  Magnified  lower  edge  of  resolution  line. 


pixel  of  the  camera  is  6.8um.  However,  when  magnification  is  taken  into  consideration, 
one  pixel  is 

- - =  0.9 107 ±0.7286%  (6.4) 

7.4664  ±0.7286% 

or 

0.9 107  um  ±  0.00663 6um .  (6.5) 

Therefore,  accuracy  of  placement  must  be  within  0.9107um.  The  most  accurate 
translation  stage  available  for  use  in  these  experiments  has  a  smallest  increment  of  0.5um. 
Straightness  of  travel  and  orientation  of  the  camera  with  respect  to  the  direction  of  travel 
of  the  stage  must  also  be  considered.  When  all  the  factors  are  combined,  it  was 
considered  that  a  test  of  mosaicking  was  impractical. 

6.3  Film  Exposure  and  Developing 

The  hologram  film  plates  used  were  8E75HD-1  Holotest  plates  manufactured  by 

Agfa.  The  specified  exposure  was  — — at  the  HeNe  laser  wavelength  to  achieve  an 

cm 

exposure  density  of  1 .  The  power  meter  used  is  a  Spectra-Physics  model  404.  It  was  not 
necessary  to  use  an  attenuator.  The  entrance  aperture  of  the  detector  head  is  specified  as 
7mm  in  diameter,  thus  the  area  is  0.385cm2.  In  order  to  increase  visibility  of  the  fringes, 
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the  weak  beam  from  the  splitter  was  used  as  the  reference  plane  wave  and  the  strong 
beam  was  used  as  the  object  wave.  Since  it  was  necessary  to  use  a  scatter  plate,  the 
object  beam  reaching  the  film  plate  was  still  weaker  them  the  reference  beam.  So,  neutral 
density  filters  were  also  used.  It  was  determined  experimentally  that  the  two  neutral 
density  filters  together  had  a  transmission  of  0.067. 

For  the  double  exposure  hologram,  the  reference  beam  was  measured  to  be 
0.006mW  and  the  object  beam  (including  the  scatter  plate)  was  measured  to  be  0.5uW. 
These  measurements  were  with  the  resolution  target  at  a  distance  of  164mm.  When  the 
neutral  density  filters  were  included,  the  reference  beam  then  became  0.402uW.  So,  the 
total  power  at  the  film  plate  was  approximately  1  .OuW.  The  exposure  time  can  then  be 
calculated  as 

15m/ 

—ft— =  5.17  5s.  (6.6) 

luW 

0.385cm2 

Since  this  was  a  double  exposure  hologram,  it  was  decided  to  reduce  the  exposure  by  a 

factor  of  ■yjl  and  so  the  exposure  then  became  4.08s.  The  actual  used  was  4.1s  because 
that  was  the  closest  setting  possible  on  the  shutter  controller.  For  the  second  exposure 
with  the  resolution  target  at  350mm,  a  similar  procedure  led  to  an  exposure  of  8.2s. 
Similarly,  for  the  depth  of  field  test  with  a  ruler  as  the  object  the  exposure  was  8.6s. 

The  developing  of  the  film  plates  was  done  as  follows.  The  developer  used  was 
Kodak  Dektol  developer  cat  #  146  4726  and  the  fixer  used  was  Kodak  fixer  cat  #  197 
1746.  The  film  plates  were  placed  in  the  developer  for  3min,  washed  in  running  water  for 
lmin,  placed  in  the  fixer  for  5min,  washed  in  running  water  for  lOmin,  and  finally  dried. 

6.4  Resolution  Determination 


The  resolution  limit  imposed  by  the  size  of  the  hologram  digitized  is  now 
addressed.  Three  different  reconstructions  are  performed  and  analyzed.  The  first  is  a 
reconstruction  from  a  hologram  made  of  the  resolution  target  at  a  distance  of  218mm 
from  the  hologram  film  plate.  The  entire  reconstruction  is  given  in  Fig.  6.4  and  Fig.  6.5 
is  a  magnified  view  of  just  the  resolution  target  itself.  In  Fig.  6.4,  the  real  and  virtual 
images  have  been  labeled.  Note  that  these  are  as  predicted  in  Fig.  2.4.  The  image  in  Fig 
6.5  has  been  rotated  so  that  the  resolution  target  matches  the  orientation  of  the  one  in  Fig. 
4.1.  The  reconstruction  is  noisy,  making  it  difficult  to  determine  the  actual  resolution. 
From  a  visual  inspection  of  Fig.  6.5,  it  appears  that  the  maximum  resolution  is  about 

Group  1  Element  4.  The  expected  resolution  was  3.38—^—  which  corresponds  to  about 


mm 


Group  1  Element  5. 


The  achievable  resolution  can  be  estimated  by  determining  the  size  of  the 
individual  pixels  that  make  up  the  reconstructed  image.  Since  the  size  of  the  lines  and 
blocks  in  the  resolution  target  are  known,  then  the  pixel  size  can  be  calculated.  The  best 
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Virtual  Image 


Real  Image 


Figure  6.4  Reconstruction  of  the  resolution  target  at  a  distance  of  218mm. 


Figure  6.5  Magnified  view  of  the  reconstructed  resolution  target 
at  a  distance  of  2 1 8mm. 
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part  of  the  resolution  target  to  use  is  the  largest  block  since  it  would  contain  the  most 
pixels  thus  reducing  the  error  in  the  result.  The  size  of  the  largest  block  is  2.5mm  X 
2.5mm.  By  zooming  in  on  this  block,  it  was  estimated  that  the  block  is  17  pixels  X  17 
pixels.  Due  to  the  noisiness  of  the  image,  an  error  estimate  of  ±  2  pixels  was  used  which 
gives  a  relative  error  of  ±  1 1.76%.  This  gives  a  pixel  size  of 

- - =  0.147 1  mm  ±  1 1.76%  .  (6.7) 

17  pixels  ±  11.76% 

Thus,  the  maximum  spatial  frequency  possible  is 

- ^ - =  3.4-^- ±11.76%  .  (6.8) 

2(0.147mm)  ±  1 1.76%  mm 

This  is  a  range  of  3.0  to  3.8-^- .  The  theoretical  value  falls  within  this  range. 


.  The  theoretical  value  falls  within  this  range. 


The  error  in  the  expected  spatial  resolution  must  also  be  evaluated.  The 
maximum  spatial  frequency  is  calculated  from  Eq.  (4.2).  The  magnification  M  and  the 
distance  d  are  the  main  factors  where  a  significant  error  can  result.  The  magnification 
and  its  error  were  given  in  Eq.  (6.3).  For  the  distance  d,  the  error  results  from  the  actual 
placement  of  the  resolution  target.  The  smallest  increment  on  the  ruler  used  in  the 
measurement  is  1mm.  However,  the  measurement  of  the  distance  was  a  little  difficult 
and  so  the  error  is  estimated  to  be  +  2mm  which  gives  a  relative  error  of  ±  0.9174%  for  a 
distance  of  164mm.  Thus,  the  calculated  resolution  and  its  error  is 

- . -  = - 1024(6.8Mm) - _  3.380-^- ±  1.646%  .  (6.9) 

2MXd  2(7.4664  ±0.7286%)(632.8wn)(218mm±  0.9174%)  mm 

In 

This  is  a  range  of  3.324  to  3.436—^— .  This  falls  well  within  the  measured  spatial 

mm 

frequency  range  given  in  Eq.  (6.8). 


The  next  resolution  test  was  at  a  distance  of  164mm.  Actually,  the  hologram  is  a 
double  exposure  hologram  of  the  resolution  target  at  164mm  and  350mm.  The  magnified 
view  of  the  reconstructed  resolution  target  is  shown  in  Fig.  6.6.  Visually,  the  maximum 


Figure  6.6  Magnified  view  of  the  reconstructed  resolution  target 
at  a  distance  of  d=  164mm. 


spatial  frequency  resolvable  is  about  Group  1  Element  6.  The  predicted  value  was 

4.493  — ,  which  corresponds  to  Group  2  Element  2. 
mm 

As  before,  the  achievable  resolution  can  be  estimated  by  determining  the  size  of 
the  individual  pixels  that  make  up  the  reconstructed  image.  Once  again,  the  largest  block 
on  the  resolution  target  is  used  to  determine  the  size  of  the  pixels.  By  zooming  in  on  the 
reconstructed  image,  it  was  estimated  that  the  largest  block  is  21  pixels  X  21  pixels. 
Again,  assuming  an  error  of  ±  2  pixels  gives  a  relative  error  of  ±  9.524%.  The  resulting 
pixel  size  is 

— 25mm  =  0.1 19mm  ±9.524%  .  (6. 10) 

21  ±9.524% 

The  corresponding  maximum  spatial  frequency  is 

- ^ - —  =  4.2 -^-±9.524%.  (6.11) 

2(0.1 1 9  mm)  ±  9524%  mm 


This  gives  a  range  of  3.8  to  4.6-^— .  The  theoretical  value  of  4.493  — falls  within  this 

mm  mm 

range. 

The  error  in  the  theoretical  value  should  be  evaluated  as  before.  The  relative  error 
of  the  distance  measurement  of  d=164±  2mm  is  ±  1.22%.  The  error  of  the  expected 
value  in  this  case  is 

_ = _ 1 024(6.8 w?Q - =  4.493— ±1.949%,  (6.12) 

2  MM  2(7.4664  ±0.7286%)(632.8nm)(164mm±  1.22%)  mm 

which  is  a  range  of  4.405  to  4.581  .  This  range  falls  within  the  range  of  the  measured 

mm 

value. 


The  final  resolution  test  is  of  the  resolution  target  at  a  distance  of  350mm.  The 
reconstruction  of  the  resolution  target  at  350mm  is  given  in  Fig  6.7.  A  visual  inspection 


Figure  6.7  Magnified  view  of  the  reconstructed  resolution  target 
at  a  distance  of  350mm. 
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of  Fig  5.10  reveals  that  the  maximum  spatial  frequency  resolvable  is  Group  0  Element  5. 

The  expected  resolution  was  2.105-^- ,  which  corresponds  to  about  Group  1  Element  1. 

mm 


Using  the  same  approach  as  before,  the  achievable  resolution  can  be  estimated  by 
determining  the  size  of  the  individual  pixels  that  make  up  the  reconstructed  image.  Again 
the  largest  block  on  the  resolution  target  is  used  to  determine  the  size  of  the  pixels.  By 
zooming  in  on  the  reconstructed  image,  it  was  estimated  that  the  largest  block  is  1 1  pixels 
X  1 1  pixels.  Again,  assuming  an  error  of  +  2  pixels,  the  relative  error  is  18.18%.  The 
resulting  pixel  size  is 

2  ^tTTTfZ  /.  i  <  n  ,  r\  rrt  ( S  1  *5  \ 

(6.13) 


11±  18.18% 

The  corresponding  maximum  spatial  frequency  is 

1  Ip 


=  0.2273mm +18.1 8% 


2(0.227 3mm)  ±18.18% 
Ip 


=  2.2 


Ip 


±18.18% 


mm 


(6.14) 


Thus,  the  range  is  1.8  to  2.6— — .  Once  again  the  theoretical  value  falls  within  the 


mm 


experimentally  determined  range. 

Following  the  same  analysis  as  previous,  the  theoretical  value  must  be  modified  to 
take  in  consideration  the  errors  in  the  parameters.  The  relative  error  in  the  distance  of 
350mm  is  +  0.5714%.  The  range  of  the  theoretical  value  is  therefore 

— C— = - 1024(6.8»m) - =  2.105-^-+  1.3%  ,  (6.15) 

2  MM  2(7.4664  ±0.7286%)(632.8«m)(350mm±  0.57 14%)  mm 

1 

which  corresponds  to  a  range  of  2.078  to  2.132 - .  This  range  falls  within  the  range  of 

mm 

the  measured  value. 


6.5  Depth  of  Field  Demonstration 

As  described  in  Chapter  4,  two  approaches  were  used  in  demonstrating  the  depth 
of  field.  The  first  approach  was  to  make  a  double  exposure  hologram  of  the  resolution 
target  at  two  different  distances.  This  is  the  same  hologram  used  in  the  previous  section 
for  the  resolution  target  at  164mm  and  350mm.  For  the  depth  of  field  demonstration,  the 
same  hologram  was  reconstructed  at  those  same  distances.  The  reconstructions  are 
shown  in  Fig.  6.8.  In  Fig.  6.8,  the  larger  resolution  target  is  the  one  located  at  164mm 
while  the  smaller  one  is  the  same  resolution  target  except  that  it  is  was  moved  to  a 
distance  of  350mm.  In  Fig  6.8a,  the  reconstruction  distance,  or  focus  distance,  was 
164mm.  In  Fig  6.8b,  the  reconstruction  distance,  or  focus  distance,  was  350mm. 
Referring  to  Fig  4.4,  notice  that  the  depth  of  field  is  much  smaller  when  the  focus  is  at 
164mm  in  comparison  with  when  the  focus  is  at  350mm.  This  effect  can  be  seen  in  Fig. 
6.8.  In  Fig.  6.8a,  the  resolution  target  at  164mm  is  in  focus  while  at  350mm  it  is  clearly 
out  of  focus.  On  the  other  hand,  in  Fig.  6.8b  the  resolution  target  at  350mm  is  in  focus, 
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but  also  the  resolution  target  at  164mm  is  still  somewhat  in  focus  although  noticeably 
degraded  compared  to  Fig.  6.8a. 


(a)  d=  164mm  (b)  d=350mm 


Figure  6.8  Images  of  the  double  exposure  hologram  reconstructed 
at  two  different  distances:  (a)  d=164mm.  (b)  d=350mm. 

The  other  demonstration  of  depth  of  field,  as  described  in  Chapter  4,  is  a 
hologram  of  a  specially  designed  “ruler”.  This  ruler  was  placed  at  an  angle  to  the 
hologram  film  plate.  The  concept  was  that  different  portions  of  the  ruler  would  come  into 
focus  depending  on  the  focus  distance.  Figure  6.9  shows  the  reconstructions  of  the  ruler 
at  distances  of  73.25mm,  106.8mm,  129.1mm,  and  185mm.  Using  Eqs.  (4.8)  and  (4.9), 
the  corresponding  depths  of  field  are  as  follows: 

d  =  73.25mm  =>  A s  =  ±5547mm 

d  =  106.8mm  =»  As  =  ±11.7 9mm  1  , 

(6.16) 

d  =  129.1mm  =>  As  =  ±  17.23mm 
d  =  1 85.0mm  =>  As  =  ±3 5.3 8mm. 

The  individual  lines  in  the  ruler  are  1.329mm  in  width.  However,  their  effective 
width  in  terms  of  distance  directly  normal  to  the  hologram  film  plate  is  1.299mm.  This 
takes  into  account  the  angle  of  the  ruler.  This  is  illustrated  in  Fig.  6.10.  By  using  the 
effective  width  of  the  lines,  the  focal  plane  and  depths  of  field  locations  can  be  found 
roughly  by  counting  the  lines. 

It  is  difficult  to  perceive  the  actual  depths  of  field  from  the  images  in  Fig.  6.9 
since  depth  of  field  is  a  subjective  measure.  However,  the  variations  in  depth  of  field  are 
definitely  noticeable.  The  focus  distance  in  Fig  6.9a  corresponds  to  the  end  of  the  ruler 
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(a)  d=73.25mm. 


(b)  d=  106. 8mm. 


Figure  6.9  Reconstruction  of  ruler  at  various  focus  distances:  (a)  d=73.25mm. 
(b)  d=  106.8mm.  (c)  d=129.1mm.  (d)  d=185.0mm. 


closest  to  the  hologram  film  plate  while  in  Fig  6.9d  the  focus  is  on  the  end  of  the  ruler 
farthest  from  the  hologram  film  plate.  As  in  the  double  exposure  hologram,  the  larger 
depth  of  field  at  longer  distances  from  the  film  plate  is  evident.  In  Fig.  6.9a,  only  a  small 
portion  of  the  ruler  is  in  focus  while  seemingly  the  entire  ruler  is  in  focus  in  Fig.  6.9d. 


Effective  width 
of  1  ruler  line  -<^ 
(to  measure  distances 
normal  to  hologram) 


Figure  6.10  Effective  width  of  ruler  line. 


6.6  Distortion  and  Noise  Considerations 

There  are  three  obvious  contributions  of  distortion  and  noise  in  the  digital 
reconstructions.  They  are  low  pass  filtering  due  to  the  pixel  size,  quantization  noise,  and 
speckle  noise.  A  description  of  each  of  these  will  now  be  given. 

A  major  contributor  of  distortion  is  the  low  pass  filtering  effect  caused  by  the 
finite  size  of  the  pixels  in  the  camera  as  discussed  in  Chapter  3.  The  case  of  the  pixel 
width  being  the  same  as  the  pixel  spacing,  as  is  true  with  the  video  camera  used  in  the 
experiments,  was  addressed.  It  was  found  that  the  result  is  a  1.96dB  reduction  in 
frequency  amplitude  for  frequencies  at  half  the  sampling  frequency. 

Quantization  noise  occurs  when  an  analog  signal  is  digitized.  The  number  of  bits 
in  the  digitization  limits  the  number  of  levels  that  are  available  to  represent  the  signal. 
This  causes  the  effect  known  as  quantization  noise.  It  is  usually  estimated  that  the  signal- 
to-noise  ratio  attributable  to  quantization  is  given  by  6dB  per  bit  [Ziemer,  et.  al.,  1983,  p. 
302].  The  digitizer  used  provides  8  bits  per  pixel;  thus  the  maximum  signal-to-noise  ratio 
is  about  48dB. 


Laser  speckle  is  another  contribution  of  noise  in  the  digital  reconstruction  of 
holograms.  Speckle  noise  occurs  due  to  variations  in  the  surface  of  the  object  and  film 
emulsion,  film  grains,  and  contaminates  in  the  environment.  It  is  this  effect  which  causes 
the  granularity  that  appears  in  reconstructed  holograms.  Characterization  of  speckles 
involves  lengthy  statistical  analysis.  The  essential  results  can  be  found  in  several  texts. 
For  instance,  in  Goodman,  1996,  p.  369,  the  speckle  size  <7S  is  given  to  be  roughly 
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where  D  is  the  size  of  the  hologram  (apparently  assuming  use  of  a  circular  area  of  the 
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hologram).  In  Vikram,  1992,  p.  62,  the  typical  speckle  size  was  given  to  be 
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where  r  is  the  radius  of  the  hologram  area  used  in  the  reconstruction.  In  Kreis,  1996,  p. 
37,  it  was  shown  that  the  average  size  of  a  speckle  for  a  uniformly  scattering  square  of 
dimensions  L  X  L  is  given  by 
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The  result  in  Eq.  (6.19)  seems  most  applicable  for  the  case  of  digital  holography  where  a 
square  area  of  the  hologram  is  typically  used  as  discussed  previously. 


The  importance  of  Eq.  (6.19)  is  that  it  suggests  noise  in  the  image  that  is  larger 
than  the  theoretical  resolution  capability  for  a  given  scanned  area  of  a  hologram  as 
derived  in  Chapter  3.  It  was  shown  [Eq.  (3.29)]  that  the  maximum  resolution  for  a 
limited  scan  area  of  a  hologram  is 


h  = 
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Thus,  the  average  speckle  size  as  given  by  Eq.  (6.19)  is  larger  than  the  theoretical 
resolution  limit  in  Eq.  (6.20)  by  a  factor  of  K . 


The  addition  of  distortion  and  noise  from  each  of  the  described  effects  will 
obviously  cause  degradation  in  the  quality  of  the  reconstructed  images.  This  will  make 
the  resolution  target  used  in  the  experiments  more  difficult  to  see  for  sets  of  lines  that  are 
close  to  the  theoretical  resolution  limit.  Perhaps  not  so  obvious,  is  the  effect  on  the  depth 
of  field.  Depth  of  field  is  determined  by  the  resolution  ability  of  a  system.  Better 
resolution  gives  a  smaller  depth  of  field.  Conversely,  worse  resolution  gives  a  larger 
depth  of  field.  This  is  clearly  noticeable  in  Fig  6.9  where  the  depths  of  field  were  larger 
than  anticipated. 


7.0  CONCLUSIONS 


7.1  Project  Results 

There  were  5  main  objectives  of  this  research  effort.  The  first  objective  was  to 
develop  the  mathematical  theory  of  digital  holography.  The  result  of  the  theory  was  that 
reconstructions  are  essentially  performed  with  the  well-known  Fourier  transform. 
Focusing  on  various  parts  of  an  object  is  achieved  by  including  in  the  transform  an 
additional  phase  term. 

The  second  objective  was  to  describe  the  requirements  and  effects  of  digitizing  a 
hologram.  Based  on  the  angles  between  the  reference  and  object  waves,  it  was  shown 
how  to  determine  the  numerical  aperture  of  the  microscope  objective  needed  to  resolve 
the  resulting  spatial  frequencies  of  the  fringes  recorded  on  the  hologram.  Additionally,  it 
was  shown  that  the  finite  size  of  the  pixels  in  the  digitizing  camera  has  a  lowpass  filtering 
effect  on  the  reconstructions.  Also,  it  was  shown  that  assembling  sub-portions  of  a 
hologram  into  one  contiguous  hologram  required  placement  accuracies  of  better  than  one 
pixel  in  order  to  achieve  any  benefit. 

The  third  objective  was  to  predict  the  spatial  resolution  possible  in 
reconstructions.  Based  on  the  mathematical  theory,  it  was  shown  how  much  of  a 
hologram  must  be  digitized  in  order  to  resolve  fine  detail  of  an  object.  The  results 
indicate  that  relatively  large  portions  of  a  hologram  must  be  digitized.  For  instance,  if  it 
is  desired  to  resolve  an  object  that  is  0.1mm  X  0.1mm  in  size  at  a  distance  of  500mm  and 
assuming  the  wavelength  is  632.8nm,  then  the  hologram  must  be  at  least  3.164mm  X 
3.164mm. 

Digitizing  an  entire  hologram  of  this  size  would  be  a  formidable  task  when 
considering  that  magnification  would  probably  be  necessary.  Flatbed  scanners  have  an 
optical  resolution  typically  around  600  dots  per  inch  (dpi).  This  is  not  sufficient  to 
resolve  the  spatial  frequencies  present  on  a  typical  hologram. 

The  amount  of  computer  hard  disk  storage  space  required  for  such  a  hologram  is 

tremendous.  If  the  angles  of  the  object  beam  and  reference  beam  are  assumed  to  be  30° 
each,  then  the  spatial  frequencies  recorded  on  the  hologram  will  be  on  the  order  of  1,580 
cycles/mm.  Assuming  sampling  at  the  Nyquist  rate,  then  100  megasamples  are  needed. 

If  the  digitization  is  8bits/sample,  then  the  necessary  storage  space  is  100  megabytes. 

Even  if  an  entire  hologram  can  be  digitized  and  stored,  a  computer  system  is  needed  that 
is  capable  of  processing  such  large  amounts  of  data  in  a  timely  fashion. 

The  mere  storage  of  the  hologram  itself  as  describe  above  is  only  part  of  the 
memory  problem.  Since  the  mathematical  operations  involved  in  the  reconstructions  are 
complex  floating-point  operations,  the  amount  of  memory  needed  to  process  the 
hologram  is  even  larger.  The  size  of  one  double-precision  complex  floating-point 
number  is  16  bytes  (8  bytes  each  for  the  real  and  imaginary  parts).  Thus,  the  hologram  in 


our  example  will  now  require  1.60  gigabytes  of  memory  to  process  the  entire  hologram  at 
one  time.  The  bottom  line  of  the  memory  issue  is  can  software  programs  process 
matrices  that  are  so  large,  is  there  enough  hard  disk  storage  and  RAM  memory  space,  is 
the  time  required  for  processing  within  tolerable  lengths,  and,  of  course,  is  such  a 
computer  system  affordable. 

The  fourth  objective  was  to  predict  the  depth  of  field  in  reconstructions.  It  was 
shown  that  the  depth  of  field  is  dependent  on  the  spatial  resolution.  It  was  seen  that  better 
resolution  implies  a  smaller  depth  of  field.  Conversely,  worse  resolution  implies  a  larger 
£  depth  of  field. 

The  final  objective  was  to  verify  the  predictions  through  experimentation.  The 
-  experiments  performed  demonstrated  close  agreement  with  the  expected  spatial 

resolution.  However,  the  depths  of  field  seen  in  the  reconstructions  were  larger  than 
anticipated.  The  primary  causes  of  this  result  are  the  distortion  and  noise  as  discussed  in 
Chapter  6  and  also  the  subjective  nature  of  perception  as  to  what  is  considered  in-focus. 

In  conclusion,  in  time,  as  the  spatial  resolution  of  optical  components  such  as 
flatbed  scanners  increases  and  the  capabilities  of  computer  systems  improves  (which  is 
currently  occurring  at  an  incredible  rate),  then  the  technique  of  digital  holography  will 
become  an  increasingly  favorable  approach  for  hologram  analysis. 

7.2  Suggestions  for  Future  Research 

There  are  several  areas  to  suggest  for  future  research.  The  mosaicking  problem  is 
deserving  of  further  attention.  The  analysis  of  this  research  effort  was  necessarily 
somewhat  simplistic,  limiting  consideration  to  only  x-y  placement  locations.  There  are  a 
number  of  papers  describing  and  measuring  the  quality  of  methods  for  mosaicking  images 
together.  These  various  methods  must  be  approached  cautiously  because  it  is  common  to 
warp  or  otherwise  process  images  to  provide  an  aesthetically  appealing  mosaic.  The 
consequences  of  such  image  processing  on  the  scientific  data  to  be  extracted  must  be 
considered. 

Y  Depth  of  field  lends  itself  to  further  research.  It  has  been  described  [Goodman, 

1996,  p.  101-107]  how  a  lens  can  be  put  in  various  locations  with  respect  to  the  hologram 
plane  during  reconstruction.  If  a  lens  can  be  simulated  and  placed  at  various  locations  in 
'  the  digital  reconstruction  process,  then  this  suggests  that  the  image  and  object  distances 

can  be  varied.  Thus,  the  depth  of  field  could  possibly  be  manipulated  to  some  degree. 

Perhaps  the  most  obvious  is  the  need  for  further  research  is  reconstructing  images 
of  three  dimensions.  Each  of  the  images  reconstructed  in  this  research  effort  was  two- 
dimensional.  Stereoscopic,  tomographic,  or  other  such  methods  could  be  applied  to 
attempt  to  extend  this  research  to  true  three-dimensional  image  reconstruction  and  data 
extraction.  Ultimately,  there  must  be  a  more  direct  computational  method  for 
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reconstructing  in  three  dimensions  rather  than  to  require  use  of  such  seemingly 
cumbersome  approaches. 


I 
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APPENDIX  A 


Mosaicking  Error  Analysis  Program 

The  program  used  to  simulate  mosaicking  for  error  analysis  is  named  encsumfp.m. 
The  following  is  a  listing  of  the  MATLAB  (by  The  MathWorks,  Inc.)  program  code: 

%  Mosaicking  Error  Analysis 

close  all 

clear  all 

whitebg('w') 

tic 

L=32768;  %  Number  of  points  in  one  segment  of  the  signal 
seg=4;  %  Number  of  segments  in  the  signal 
f=l/15;  %  Frequency  of  the  signal 

iter=500;  %  The  number  of  iterations  (for  statistical  'purposes) 
maxdev=2.5*pi;  %  This  is  the  maximum  phase  angle  deviation  that  will  be  tested 
phasenum=5;  %  The  number  of  steps  leading  up  to  the  max  phase  angle  deviation 
points=  1048576;  %  The  signal  is  zero  padded  to  this  number  of  points 
maxx=30;  %  The  encircled  power  is  computed  up  to  this  radius  (in  pixels) 

%  This  set  of  lines  converts  the  variables  from  numbers  to  string  for  use  in  labelling 
graphs 

Lstr=num2str(L) ; 
segstr=num2str(seg); 
iterstr=num2str(iter) ; 
phasenumstr=num2str(phasenum+ 1 ) ; 
fstr=num2str(f); 
pointsstr=num2str(points,8) ; 
maxdevstr=num2str(maxdev); 

%  These  lines  create  the  first  segment  of  the  signal 
i=l:L; 

yseg0=cos(2*pi*f*i); 

%  These  lines  perform  a  Fourier  transform,  discard  half  of  the  result,  and  then  square  it 
freq=fft(ysegO, points) ; 
freq=freq(  1  :round(points/2)); 
freq=(abs(freq)).A2; 

%  These  lines  locate  the  maximum  of  the  energy  computed  above,  determines  the 

encircled  energy  of  increasingly  larger  radii,  and  plots  the  results 

[z,c]=max(freq); 

energy=sum(freq) ; 

en=zeros(  1  ,maxx+ 1 ) ; 

for  p=0:maxx 

en(p+ 1  )=en(p+ 1  )+sum(freq((c-p) :  (c+p))) ; 

end 

prcnten  1  seg=en/energy ; 
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plot(0:maxx,prcnten  1  seg,'k— ’) 

for  n=0:phasenum  %  This  for  loop  increments  the  maximum  phase  deviation 
n 

%  These  lines  initialize  variables  to  0 
totalphi=0; 
en=zeros(  1  ,maxx+l); 
freq=zeros(  1  ,round(points/2)); 
for  j=l:iter  %  This  for  loop  iterates  the  simulation 
segphi=zeros(l,seg);  %  Initialize  segphi  to  0 

%  This  set  of  lines  creates  the  signal  with  each  segment  given  a  random  \ 

phase  with  respect  to  the  first  segment 
y=ysegO; 
for  m=2:seg 

segphi(m)=segphi(m-l)+(rand(l)-0.5)*n/phasenum*maxdev;  % 

This  computes  a  random  phase  for  each  segment  adds  it  to  the 
phase  of  the  previous  segment 
yseg=cos(2*pi  *f*  (i+(m- 1  )*L)+segphi(m)); 
y=[y,yseg]; 
end 

%  These  lines  perform  a  Fourier  transform,  discard  half  of  the  result,  and 

then  square  it 

newfreq=fft(y,points) ; 

newfreq=newfreq(  1  :round(points/2)); 

newfreq=(abs(newfreq)).A2; 

freq=freq+newffeq;  %  The  result  from  the  above  lines  is  summed  for  each 
iteration  and  averaged  later 

totalphi=totalphi+sum(segphi);  %  The  phases  are  summed  for  each 
iteration  and  averaged  later.  The  average  should  be  about  zero, 
end 

freq=freq/iter; 

%  These  lines  compute  the  encircled  energy  for  increasingly  larger  radii 
if  n==0 

[z,c]=max(freq);  %  This  determines  the  location  of  the  maximum  energy 
value  for  the  ideal  signal  and  uses  it  in  all  subsequent  calculations  of  the 
encircled  energy 
end 

energy=sum(fr eq) ; 
en=zeros(  1  ,maxx+ 1 ); 
for  p=0:maxx 

en(p+ 1  )=en(p+ 1  )+sum(ff eq((c-p):  (c+p))) ; 
end 

prcnten=en/energy; 

%  The  rest  of  the  lines  plot  the  results  and  put  titles  and  labels  on  the  graphs 

figure(l) 

hold  on 
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grid  on 
zoom  on 

xlabel('Distance  from  Center  Frequency  (pixels)') 
ylabel('Percent  of  Total  Energy') 

strtitle=['Seg  L='  Lstr  #  of  seg='  segstr  f='  fstr ',  iter='  iterstr ',  maxdev=' 

maxdevstr ',  0  pad  to  '  pointsstr]; 

title(strtitle) 

plot(0:maxx,prcnten,'k') 
hold  off 

[  figure(n+2) 

plot(0:maxx,prcnten  1  seg,'k— ') 
hold  on 

*  zoom  on 

plot(0:maxx,prcnten,'k') 
hold  off 

xlabel('Distance  from  Center  Frequency  (pixels)') 
ylabel('Percent  of  Total  Energy') 
phidevstr=num2str(n/phasenum*maxdev); 
t=toc 

tstr=num2str(t); 

strtitle=['Seg  L='  Lstr ',  #  of  seg=’  segstr ',  f='  fstr ',  iter='  iterstr ',  phidev=' 

phidevstr 0  pad  to  '  pointsstr]; 

title(strtitle) 

avgphi(n+ 1  )=totalphi/iter/seg; 
end 
t=toc 

tstr=num2str(t); 

x=maxdev*(0:phasenum)/phasenum; 

figure(n+3) 

zoom  on 

plot(x,avgphi,'k') 

xlabel('Maximium  Deviation  from  0  (radians)') 
ylabel('Average  Phi  (radians)’) 

*  figure(n+4) 
zoom  on 
plot(y,'k') 
figure(n+5) 
zoom  on 
plot(freq,'k') 
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APPENDIX  B 


Mosaicking  Program  and  Subprograms 

The  following  source  code  was  written  in  IDL  (a  product  of  Research  Systems, 
Inc).  The  source  code  is  not  extensively  commented.  Its  inclusion  is  intended  mainly  for 
reference  and  documentation.  Some  of  the  source  code  was  obtained  while  at  a  training 
class  given  by  Research  Systems,  Inc.  During  programming,  it  was  decided  that  writing 
one  large  program  was  not  the  best  approach  because  there  would  be  a  desire  to 
experiment  with  variations  in  the  source  code.  Therefore,  different  functions  were  broken 
out  into  individual  programs  so  that  when  changes  were  made,  only  one  program  had  to 
be  modified.  The  programs  included  in  this  appendix  are  the  following: 

mosaic.pro 

error.pro 

getseam.pro 

mosdims.pro 

picklist.pro 

read_mos.pro 

seaminfo.pro 

The  source  code  for  mosaic.pro  is  the  following: 

.  **** ************************************************** *********** 

9 

. ********* ************************************************* ******* 

9 

pro  zoom_event,  event 

Widget_Control,  event.top,  get_uvalue=info,  /No_Copy 
widget_control,  info.imageid,  get_uvalue=images,  /no_copy 
possibleEventTypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  NE  'DOWN'  THEN  RETURN 

info.mousex=event.x 

info.mousey=event.y 

;Change  the  event  handler  for  the  draw  widget  and  turn  MOTION  events  ON. 
Widget_Control,  event. id,  Event_Pro=’zoommove_event’,  Draw_Motion_Events=l 

;  Take  over  color  index  1  for  the  zoom  box  drawing  color.  Store  the  current  (r,g,b)  values 

;  for  color  index  1  so  you  can  restore  the  current  colors  after  the  zoom  box  is  drawn. 

TVLct,  r,  g,  b,  /Get 

info.r_old  =  r(l) 

info.g_old  =  g(l) 

info.b_old  =  b(l) 


widget_control,  info.imageid,  set_uvalue=images,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end ;  of  zoom_event 
> 

pro  zoommove_event, event 

;  This  event  handler  continuously  draws  and  erases  the  image  box  until  it  receives  an  UP 
;  event  from  the  draw  widget.  Then  it  turns  draw  widget  motion  events  OFF  and  changes 
;  the  event  handler  for  the  draw  widget  back  to  zoom_PROCESS_EVENTS. 

;  Get  the  info  structure  out  of  the  top-level  base. 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 

widget_control,  info.imageid,  get_uvalue=images,  /no_copy 

;  What  type  of  an  event  is  this? 

possibleEventT ypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  EQ  'UP'  THEN  BEGIN 

;  If  this  is  an  UP  event,  you  need  to  erase  the  zoombox,  restore  the  user's  color 
;  table,  turn  motion  events  OFF,  set  the  draw  widget’s  event  handler  back  to 
;  zoom_PROCESS_E VENTS,  and  draw  the  "zoomed"  plot  in  both  the  draw 
;  widget  and  the  pixmap. 

info.zoomxlow=info.mousex 

info.zoomxhigh=event.x 

info.zoomylow=info.mousey 

info.zoomyhigh=event.y 

;  Make  sure  the  x  and  y  values  are  ordered  correctly 
IF  info.zoomxlow  GT  info.zoomxhigh  THEN  begin 
info.zoomxlow=event.x 
info.zoomxhigh=info.mousex 

endif 

IF  info.zoomylow  GT  info.zoomyhigh  THEN  begin 
info.zoomylow=event.y 
info.zoomyhigh=info.mousey 


endif 


;  Scale  the  coordinates 

info.zoomxlow=round(info.zoomxlow*info.xscale) 

info.zoomxhigh=round(info.zoomxhigh*info.xscale) 

info.zoomylow=round(info.zoomylow*info.yscale) 

info.zoomyhigh=round(info.zoomyhigh*info.yscale) 

WSet,  info.widl 

Device,  Copy  =  [0, 0, 512, 512, 0, 0,  info.boxwid] 

;  Restore  the  user's  color  table 

TVLct,  info.r_old,  info.g_old,  info.b_old,  1 

;  Turn  motion  events  and  button  events  off 

Widget_Control,  event.id,  Draw_Motion_Events=0,  draw_button_events=0 
display_zoom,  info,  images 

widget_control,  info.imageid,  set_uvalue=images,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location  and  return. 
Widget_Control,  event. top,  Set_UValue=info,  /No_Copy 
RETURN 

ENDIF  ;  thisEvent  =  UP 

;  Most  of  the  action  in  this  event  handler  occurs  here  while  we  are  waiting  for  an  UP 
;  event  to  occur.  As  long  as  we  don't  get  it,  keep  erasing  the  old  zoom  box  and  drawing 
;  new  one. 

;  Erase  the  old  zoom  box. 

WSet,  info.widl 

Device,  Copy  =  [0, 0,  512,  512,  0,  0,  info.boxwid] 

;  Update  the  dynamic  comer  of  the  zoom  box  to  the  current  cursor  location, 
newx  =  event.x 
newy  =  event.y 

;  Load  a  color  in  color  index  1  to  draw  the  zoom  box  with. 

TVLct,  255B,  0B,  0B,  1 

;  Draw  the  zoom  box. 

PlotS,  [info.mousex,  info.mousex,  newx,  newx,  info.mousex],  $ 

[info.mousey,  newy,  newy,  info.mousey,  info.mousey],  /Device,  Color=l 

widget_control,  info.imageid,  set_uvalue=images,  /no_copy 


;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end  ;  of  zoommove_event 

.  ***************************************************************** 

5 

pro  select_event,  event 

Widget_Control,  event.top,  get_uvalue=info,  /No_Copy 
widget_control,  info.imageid,  get_uvalue=images,  /no_copy 
possibleEventTypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  NE  'DOWN'  THEN  RETURN 

;  Obtain  the  first  mouse  coordinates 

info.mousex=event.x 

info.mousey=event.y 

;  Determine  the  selected  image  from  the  mouse  coordinates  if  the  image  was  not 
;  selected  from  the  List, 
if  info.list  eq  0  then  begin 

for  i=0,info.imnum-l  do  begin 

if  (info.mousex*info.xscale  ge  images(i).x)  and  $ 

(info.mousex*info.xscale  le  (images(i).x+info.subxsize))  and  $ 
(info.mousey*info.yscale  ge  images(i).y)  and  $ 
(info.mousey*info.yscale  le  (images(i).y+info.subysize))  then  $ 
info.imsel=i ;  The  HIGHEST  image  #  in  the  area  will  be  selected 

endfor 

endif 

;  Take  over  color  index  1  for  the  zoom  box  drawing  color.  Store  the  current  (r,g,b)  values 

;  for  color  index  1  so  you  can  restore  the  current  colors  after  the  zoom  box  is  drawn. 

TVLct,  r,  g,  b,  /Get 

info.r_old  =  r(l) 

info.g_old  =  g(l) 

info.b_old  =  b(l) 

;  Load  a  color  in  color  index  1  to  draw  the  zoom  box  with. 

TVLct,  255B,  OB,  OB,  1 

;  Change  the  event  handler  for  the  draw  widget  and  turn  MOTION  events  ON. 
Widget_Control,  event. id,  Event_Pro='move_event',  Draw_Motion_Events=l 

widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
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;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end ;  of  select_event 

•  5jc  5jC  jjl  5|C  ^  5{C  jjC  /{»  jjC  jjv  ){*  5jC  !^C  ?jC  5j>  ^  ^  2-fC  5jC  ijC  JjC  5jC  5jC  *fC  ?fC  *$%  ^  ^j>  5j>  ^jc  ^  ?jC  5f%  ^  ^ 

9 

pro  move_event, event 

;  This  event  handler  continuously  draws  and  erases  the  image  box  until  it  receives  an  UP 
;  event  from  the  draw  widget.  Then  it  turns  draw  widget  motion  events  OFF  and  changes 
;  the  event  handler  for  the  draw  widget  back  to  zoom_PROCESS_EVENTS. 

;  Get  the  info  structure  out  of  the  top-level  base. 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 

widget_control,  info.imageid,  get_uvalue=images,  /no_copy 

;  What  type  of  an  event  is  this? 

possibleEventTypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL*  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  EQ  'UP'  THEN  BEGIN 

;  If  this  is  an  UP  event,  you  need  to  erase  the  zoombox,  restore  the  user's  color 
;  table,  turn  motion  events  OFF,  set  the  draw  widget's  event  handler  back  to 
;  zoom_PROCESS_EVENTS,  and  draw  the  "zoomed"  plot  in  both  the  draw 
;  widget  and  the  pixmap. 

;  store  the  new  coordinates  for  the  selected  image 
delx  =  (event.x-info.mousex)*info.xscale 
dely  =  (event.y-info.mousey)*info.yscale 
images(info.imsel).x=images(info.imsel).x+delx 
images(info.imsel).y=images(info.imsel).y+dely 

;  Erase  the  zoombox  one  final  time  by  copying  the  plot  from  the  pixmap. 

WSet,  info.widl 

Device,  Copy  =  [0, 0,  512,  512,  0,  0,  info.boxwid] 

;  Restore  the  user's  color  table. 

TVLct,  info.r_old,  info.g_old,  info.b_old,  1 

;  Turn  motion  events  and  button  events  off 

Widget_Control,  event.id,  Draw_Motion_Events=0,  draw_button_events=0 


display_mosaic,  info,  images 

widget_control,  info.imageid,  set_uvalue=images,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location  and  return. 
Widget_Control,  event. top,  Set_UValue=info,  /No_Copy 
RETURN 

ENDEF  ;  thisEvent  =  UP 

;  Most  of  the  action  in  this  event  handler  occurs  here  while  we  are  waiting  for  an  UP 
;  event  to  occur.  As  long  as  we  don't  get  it,  keep  erasing  the  old  zoom  box  and  drawing 
;  new  one. 

;  Erase  the  old  zoom  box. 

WSet,  info.widl 

Device,  Copy  =  [0,  0,  512,  512,  0,  0,  info.boxwid] 

;  Update  the  dynamic  comer  of  the  zoom  box  to  the  current  cursor  location, 
delx  =  event. x-info.mousex 
dely  =  event.y-info.mousey 

;  Draw  the  zoom  box. 

PlotS,  [round(images(info.imsel).x/info.xscale+delx),  $ 

round((images(info.imsel).x+info.subxsize)/info.xscale+delx),  $ 
round((images(info.imsel).x+info.subxsize)/info.xscale+delx),  $ 
round(images(info.imsel).x/info.xscale+delx),  $ 
round(images(info.imsel).x/info.xscale+delx)],  $ 
[round(images(info.imsel).y/info.yscale+dely),  $ 
round(images(info.imsel).y/info.yscale+dely),  $ 
round((images(info.imsel).y+info.subysize)/info.yscale+dely),  $ 
round((images(info.imsel).y+info.subysize)/info.yscale+dely),  $ 
round(images(info.imsel).y/info.yscale+dely)],  /Device,  Color=l 

widget_control,  info.imageid,  set_uvalue=images,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end  ;  of  move_event 

. ***************************************************************** 

5 

;  Display  Mosaic  procedure 
pro  display_mosaic,  info,  images 


mosaic=bytarr(info.xmos,info.ymos) 


for  i=0,info.imnum-l  do  begin 

mosaic(images(i).x,images(i).y)=images(i).img 

endfor 

wset,  info.widl 
erase 

tvscl,  congrid(mosaic,512,512) 

widget_control,info.mosaicid,set_uvalue=mosaic,  /no_copy 
end  ;  of  Display  Mosaic 

. *************** **************************************** ********** 

5 

f 

;  Display  Zoom  procedure 
pro  display_zoom,  info,  images 

mosaic=bytarr(info.xmos,info.ymos) 
for  i=0,info.imnum-l  do  begin 

mosaic(images(i).x,images(i).y)=images(i).img 

endfor 

wset,  info.widl 
erase 

tvscl, congrid(mosaic(info.zoomxlow:info.zoomxhigh,info.zoomylow:info.zoomyhigh),  $ 

512,512) 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
end  ;  of  Display  Zoom 

.  *:£*************************************************************** 

PRO  mosaic_Event,  event 

widget_control,event.top,get_uvalue=info,  /No_Copy 
widget_control,event.id,get_uvalue=value 
widget_control,info.mosaicid,get_uvalue=mosaic,  /No_Copy 
widget_control,info.imageid,get_uvalue=images,  /No_Copy 

CASE  value  OF 

'MOSDIM'  :  begin 

widget_control,event.top,tlb_get_offset=offsets  ;obtain  mosaic  X  &  Y  size 
stuff=mosdims(offsets(0)+50,offsets(l)+75,  info.xmos,  info.ymos,  $ 
info.subxsize,  info.subysize,info.subnum) 
if  stuff.xsize  ne  0  then  begin 
info.xmos=stuff.xsize 
info  .ymos=stuff .ysize 
info.subnum=stuff.subnum 
info.subxsize=stuff.subxsize 
info.subysize=stuff.subysize 
mosaic=bytarr(info.xmos,  info.ymos) 

images={img:bytarr(info.subxsize,info.subysize),x:0,y:0,name:"} 
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images=replicate(images,info.subnum) 
info.xscale=info.xmos/5 12. 
info.yscale=info.ymos/5 12. 

Window,  /Free,  XSize=512,  YSize=512,  /Pixmap 
info.boxwid  =  ID.Window 

widget_control,  info.mosaicid,  set_uvalue=mosaic,  /no_copy 
widget_control,info.imageid,set_uvalue=images,  /no_copy 
widget_control,event.top,set_uvalue=info,  /No_Copy 

endif 

end 

'LOAD'  :begin 

widget_control,event.top,tlb_get_offset=offsets  ;load  images 
pick=dialog_pickfile(/read,  path=info.impath,  get_path=impath) 
info.impath=impath 

;determine  the  image  type 
pick=strlowcase(pick) 
type=strmid(pick,strlen(pick)-3,3) 
if  type  eq  'byt'  then  imtype=0 
if  type  eq  'tif  then  imtype=l 
if  type  eq  'byt'  or  type  eq  'tif  then  begin 
images(info.imnum).name=pick 
info.imtype=imtype 
imnum=info.imnum 

images(info.imnum).img=read_mos(images(info.imnum).name, 
info.imtype,  info.subxsize,  info.subysize,  info.widl) 
info.imnum=info.imnum+l 
display_mosaic,  info,  images 

widget_control, info.mosaicid, set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
widget_control,  event.top,  set_uvalue=info,  /No_Copy 
endif  else  nothing=error(offsets(0)+50,offsets(l)+75) 

end 

'SELECT'  :begin 

widget_control,  info.imageid,  draw_Button_Events=l,  $ 
event_pro='select_event’ 
info.list=0 
wset,info.boxwid 

device,  copy  =  [0, 0,  512,  512,  0,  0,  info.widl] 
widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 


widget_control,event.top,tlb_get_offset=offsets  ;obtain  new  coordinates 
info.list=l 

names=strarr(info.subnum) 
for  i=0,info.subnum-l  do  begin 
names(i)=images(i).name 

endfor 

pick=picklist(offsets(0)+50,offsets(l)+75,info.subnum,info.imsel, names) 
info.imsel=pick 

widget_control,  info.imageid,  draw_Button_Events=l,  $ 
event_pro='select_event' 

;  Make  a  copy  of  the  current  display 
wset,info.boxwid 

device,  copy  =  [0, 0,  512,  512, 0, 0,  info.widl] 

;  Take  over  color  index  1  for  the  zoom  box  drawing  color.  Store  the 
;  current  (r,g,b)  values  for  color  index  1  so  you  can  restore  the  current 
;  colors  after  the  zoom  box  is  drawn. 

TVLct,  r,  g,  b,  /Get 
info.r_old  =  r(l) 
info.g_old  =  g(l) 
info.b_old  =  b(l) 

;  Load  a  color  in  color  index  1  to  draw  the  zoom  box  with. 

TVLct,  255B,  OB,  OB,  1 

;  Set  the  current  window  to  the  display, 
wset, info.widl 

;  Draw  the  box. 

PlotS,  [round(images(info.imsel).x/info.xscale),  $ 

round((images(info.imsel).x+info.subxsize)/info.xscale),  $ 
round((images(info.imsel).x+info.subxsize)/info.xscale),  $ 
round(images(info.imsel).x/info.xscale),  $ 
round(images(info.imsel).x/info.xscale)] ,  $ 
[round(images(info.imsel).y/info.yscale),  $ 
round(images(info.imsel).y/info.yscale),  $ 
round((images(info.imsel).y+info.subysize)/info.yscale),  $ 
round((images(info.imsel).y+info.subysize)/info.yscale),  $ 
round(images(info.imsel).y/info.yscale)],  /Device,  Color=l 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event. top,  Set_UValue=info,  /No_Copy 
end 


'SHOW' : begin 

mosaic=mosaic*0. 

mosaic(images(info.imsel).x,images(info.imsel).y)=  $ 
images(info.imsel)  .img 
wset,  info.widl 
erase 

tvscl,  congrid(mosaic,512,512) 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

'SEAM'  :begin 

widget_control, event.top, tlb_get_offset=offsets  ;obtain  new  coordinates 

info.seam=getseam(offsets(0)+50,offsets(  1  )+75) 

nothing=seaminfo(offsets(0)+50,  offsets(  1  )+75) 

cursor,  corx,  cory,  /down,  /device 

info.corx=round(corx  *info .  xscale) 

info.cory=round(cory*info.yscale) 

print,'corx=',info.corx,'cory=',info.cory 

widget_control,info.mosaicid,set_uvalue=mosaic, /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

'CORMATCH'  :begin 

widget_control , event  .top  ,tlb_get_offset=offsets 

cormax=0 

xloc=0 

yloc=0 

;  corx  &  cory  are  the  user  selected  cursor  points.  The  user  selects  a  point 
;  on  the  mosaic  for  matching.  The  user  selects  whether  to  use  a  horizontal 
;  or  vertical  line  for  matching.  The  length  of  the  line  is  defined  by  the 
;  variable  length.  The  result  is  stored  in  corvee  1.  The  2  for  loops  take  this 
;  line  and  looks  for  the  best  match  within  the  selected  image  defined  by 
;  info.imsel.  The  variable  range  defines  the  number  of  lines  (either 
;  horizontal  or  vertical)  to  search  through  within  the  selected  image.  The 
;  indices  i  and  j  define  the  coordinates  within  the  selected  image  to  search. 

;  Equally  sized  lines  at  these  coordinates  are  iteratively  stored  in  corvec2 
;  and  compared  with  corvee  1.  xloc  &  yloc  are  the  coordinates  of  the  best 
;  match. 

length=256;  line  length  to  use  in  matching 
range=100 

if  info.seam  eq  0  then  begin;  horizontal  matching 

corvecl=mosaic(info.corx:(info.corx+(length-l)),info.cory) 
for  i=(info.subxsize/2-range),(info.subxsize/2+range)  do  begin 


for  j=0,100  do  begin;  #  of  lines  to  search 

corvec2=images(info.imsel)  .img(i :  (i+(length- 1 ))  ,j ) 
coramount=correlate(corvec  1  ,corvec2) 
if  coramount  gt  cormax  then  begin 
cormax=coramount 
xloc=i 
yloc=j 

endif 

endfor 

endfor  f 

endif  else  begin;  vertical  matching 

corvecl=mosaic(info.corx,info.cory:(info.cory+(length-l))) 
for  i=0, 100  do  begin;  #  of  lines  to  search  * 

for  j=(info.subysize/2-range),(info.subysize/2+range)  $ 
do  begin 

corvec2=images(info.imsel)  .img(i  ,j :  (j +(length- 1 ))) 
coramount=correlate(corvec  1  ,corvec2) 
if  coramount  gt  cormax  then  begin 
cormax=coramount 
xloc=i 
yloc=j 

endif 

endfor 

endfor 

endelse 

print, 'xloc=',xloc,,yloc=',yloc 
if  (info.corx-xloc)  It  0  or  (info.cory-yloc)  It  0  $ 

or  (info.corx+(info.subxsize-xloc))  gt  info.xmos  $ 
or  (info.cory+(info.subysize-yloc))  gt  info.ymos  then  begin 
nothing=error(offsets(0)+50,offsets(  1  )+75) 
endif  else  begin 

;  moves  selected  image  to  the  matching  location 

images(info.imsel).x=info.corx-xloc 

images(info.imsel).y=info.cory-yloc 

endelse 

mosaic=bytarr(info.xmos, info.ymos) 
for  i=0,info.imnum-l  do  begin 

mosaic(images(i).x,images(i).y)=images(i).img 

endfor 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event. top,  Set_UValue=info,  /No_Copy 
end 

'DISPLAY' :begin 
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display_mosaic,  info,  images 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

'ZOOM'  :begin 

widget_control,  info.imageid,  draw_Button_Events=l,  $ 
event_pro='zoom_event' 
wset.info.boxwid 

A  device,  copy  =  [0, 0,  5 1 2,  5 1 2,  0,  0,  info. wid  1  ] 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
■*  Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end 

'DISPZOOM' : begin 

display_zoom,  info,  images 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

’SAVE’  :begin 

mosaic=reverse(mosaic,2) 

pick=dia!og_pickfile(/write,  path=info.impath,  get_path=impath) 
info .  impath=impath 
tiff_write, pick, mosaic 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

’SAVEDAT’  :begin 

stuff=dialog_pickfile(/write,  path='c:Y,  filter-*. dat') 
save,  mosaic,  images,  info,  filename=stuff 
widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
*  widget_control,  info.imageid,  set_uvalue=images,  /no_copy 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

'  ’RESTORE’  :begin 

widl=info.widl 

stuff=dialog_pickfile(/read,  path='c:Y,  filter='*.dat') 

restore,  stuff 

info.widl=widl 

Window,  /Free,  XSize=512,  YSize=512,  /Pixmap 
info.boxwid  =  ID.Window 
display_mosaic,  info,  images 

widget_control,info.mosaicid,set_uvalue=mosaic,  /No_Copy 
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widget_control,  info.imageid,  set_uvalue=images,  /no_copy 
Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
end 

'EXIT'  :Widget_Control,  event.top,  /Destroy 
ENDCASE 

END ;  of  mosaic_Event 

.  ***************************************************************** 

? 

PRO  mosaic 

;  images  are  transferred  from  bottom  to  top,  i.e.  the  row  with  a  0  subscript  is  drawn  on  the 

;  bottom 

!order=0 

device,  retain=2 

base  =  Widget JBase(Title='Image  Mosaicking’,  Column=2) 
drawbase=widget_base(base,column=2) 
btnbase=widget_base(base,column=  1 ) 

btn  l=widget_button(btnbase,value= 'Define  Mosaic  Dimensions',uvalue='MOSDIM') 
btn2=widget_button(btnbase,value='Load  Image',uvalue=’LOAD') 
btn3=widget_button(btnbase,value='Select  Image  with  Mouse',uvalue='SELECT') 
btn4=widget_button(btnbase,value='Select  Image  from  List',uvalue='LIST') 
btn5=widget_button(btnbase,value='Show  Current  Image', uvalue='SHOW') 
btnl  l=widget_button(btnbase,value='Define  Seam',uvalue='SEAM') 
btn  1 2=widget_button(btnbase,value=’Match  Image', uvalue='CORMATCH') 
btn6=widget_button(btnbase, value= 'Display  Mosaic’,uvalue='DISPLAY') 
btn7=widget_button(btnbase,value='Zoom',uvalue='ZOOM') 

btnlO=widget_button(btnbase,value='Display  Mosaic  Zoom  Area',uvalue='DISPZOOM') 

btn8=widget_button(btnbase,value='Save  Mosaic  (TIF)',uvalue='SAVE') 

btn  1 3=widget_button(btnbase,value=’Save  Progress',uvalue='S  AVEDAT') 

btnl4=widget_button(btnbase,value='Restore  Progess',uvalue='RESTORE') 

btn9=widget_button(btnbase,value='Exit',uvalue='EXIT') 

draw  1  =widget_draw(drawbase,xsize=5 1 2,ysize=5 1 2) 

Widget_Control,  base,  /Realize 
widget_control,draw  1  ,get_value=wid  1 

info={xmos:2048,ymos:2048,zoomxhigh:50,zoomxlow:40,zoomyhigh:50,  $ 

zoomylow:40,widl  rwidl  ,mosaicid:drawbase,imtype:0,imageid;draw  1 ,  $ 
disptype:  1  ,imnum:0,imsel:0,subnum:4,subxsize:  1024,subysize:  1024,seam:0,  $ 
corx:0,cory:0,r_old:255b,g_old:255b,b_old:255b,mousex:0,mousey:0,  $ 
boxwid:0b,xscale:1.0,yscale:1.0,list:0,impath:'c:\'} 
mosaic=Ob 
images=0b 
loadct,0 
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wset,widl 

erase 

widget_control,base,set_uvalue=info,  /No_Copy 
widget_control,drawbase,set_uvalue=mosaic,  /No_Copy 
widget_control, drawl  ,set_uvalue=images,  /No_Copy 
XManager,  'mosaic',  base,  Event='mosaic_Event' 

END ;  of  mosaic 

. H^He*:*:**#********************************************************* 
y 

.  ***************************************************************** 


The  source  code  for  err  or. pro  is  the  following: 

.  ***************************************************************** 
y 

. ***************************************************************** 
y 

Pro  error_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET_BUTTON'  THEN  BEGIN 
Widget_Control,  event.top,  /Destroy 

endif 

End  ;  of  error_Event 

.  ***************************************************************** 


Function  error,  xoff,  yoff 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Program  Error!',  xsize=150) 
ok  =  Widget_Button(base,  Value='OK') 

Widget_Control,  base,  /Realize 
XManager,  'error',  base,  /Modal 
return, 0 
End ;  of  error 

.  ** ********************************************************* ****** 
y 

.  He**** ***************************************************** ******* 


The  source  code  for  gets  earn,  pro  is  the  following: 

.  ***************************************************************** 
. ***************************************************************** 
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Pro  getseam_Event,  event 


name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET_BUTTON'  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

;  O=horizontal,  l=vertical 
’HORIZONTAL'  :  Begin 
data  =  0 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

'VERTICAL'  :  Begin 
data  =  1 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

ENDCASE 


ENDIF 

End  ;  of  getseam_Event 
Function  getseam,  xoff,  yoff 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 
Title='Click  on  Seam  Direction') 
horiz  =  Widget_Button(base,  Value= 'HORIZONTAL') 
vert  =  Widget_Button(base,  Value='VERTICAL') 
Widget_Control,  base,  /Realize 
handle  =  handle_create() 
info=  { handle:  handle } 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'getseam',  base,  /Modal 
Handle_Value,  handle,  data,  /No_Copy 
Handle_Free,  handle 
Return,  data 
End ;  of  getseam 

9 
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The  source  code  for  mosdims.pro  is  the  following: 

.  *** ** ***************************************************** ******* 
5 

.  ***************************************************************** 


Pro  mosdims_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  ’WIDGET_BUTTON'  THEN  BEGIN 
h  Widget_Control,  event,  top,  Get_UValue=info,  /No_Copy 

Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

'CANCEL'  :  Begin 
data={xsize:0} 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

'ACCEPT'  :  Begin 

Widget_Control,  info.xsize_id,  Get_Value=xsize 
Widget_Control,  info.ysize_id,  Get_Value=ysize 
Widget_Control,  info.subxsize_id,  Get_Value=subxsize 
Widget_Control,  info.subysize_id,  Get_Value=subysize 
Widget_Control,  info.subnum_id,  Get_Value=subnum 
data  =  {xsizeixsize,  ysize:ysize,  subnum:subnum,  $ 
subxsizeisubxsize,  subysizeisubysize} 
Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
End 

ENDCASE 


ENDIF 

END  ;  mosdims_Event 

^  .sic**************************************************************** 

* 

Function  mosdims,  xoff,  yoff,  currentx,  currenty,  currentsubx,  currentsuby,  $ 
currentsubnum 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Enter  Mosaic  Information') 
mosbase  =  Widget_Base(base,  Column=l) 
filelabel  =  Widget_Label(mosbase,  Value='Input  Mosaic  Size') 
xsize_id  =  CW_Field(mosbase,  Title=’X  Size: ',  Value=currentx,  $ 

/Integer,  /Retum_Events) 

ysize_id  =  CW_Field(mosbase,  Title='Y  Size: ',  Value=currenty,  $ 
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/Integer,  /Return_Events) 
subnumbase=Widget_Base(base,  Column=l) 

fllelabel2  =  Widget_Label(subnumbase,  Value=’Input  #  of  Sub-Images') 

subnum_id  =  CW_Field(subnumbase,  Title='Number:’,  Value=currentsubnum,  /Integer,  $ 

/Retum_Events) 

subsizebase=Widget_Base(base,  Column=  1 ) 

filelabeB  =  Widget_Label(subsizebase,  Value='Input  Sub-Image  Size') 

subxsize_id  =  CW_Field(subsizebase,  Title='X  Size: ',  Value=currentsubx,  /Integer,  $ 

/Retum_Events) 

subysizejd  =  CW_Field(subsizebase,  Title='Y  Size:  ’,  Value=currentsuby,  /Integer,  $  4 

/Retum_Events) 

btnbase  =  Widget_Base(base,  Row=l,  /Frame) 
cancel  =  Widget_Button(btnbase,  Value='CANCEL’) 
done  =  Widget_Button(btnbase,  Value=’ ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {xsize_id:xsize_id,  ysize_id:ysize_id,  handle: handle,  currentx:currentx,  $ 
currentyxurrenty,  subnum_id:subnum_id,  $  subxsize_id:subxsize_id,  $ 
subysize_id :  suby  size_id } 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'mosdims',  base,  /Modal 

;Must  use  handle  because  all  widgets  are  destroyed  after  Xmanager. 

Handle_Value,  handle,  data,  /No_Copy 

Handle_Free,  handle 
Return,  data 
End ;  of  mosdims 

.  * ******************************************************* ********* 
.***************************************************************** 


The  source  code  for  picklist.pro  is  the  following: 

.  * ** ****************************************************** ******** 
.  ***************************************************************** 

Pro  picklist_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  ’WIDGET.LIST’  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
data=event.index 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
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Widget_Control,  event.top,  set_UValue=info,  /No_Copy 
Widget_Control,  event.top,  /Destroy 

ENDIF 

END ;  of  picklist_Event 

.  *************************** ******************** ****************** 

9 

Function  picklist,  xoff,  yoff,  number,  imsel,  names 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  Title='Image  List') 
list=widget_list(base,  value=names,  ysize=number) 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 
info=  { handle:handle,listid:list } 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
widget_control,  list,  set_uvalue=names,  /no_copy 
XManager,  'picklist',  base,  /Modal 

;Must  use  handle  because  all  widgets  are  destroyed  after  Xmanager. 

HandleJValue,  handle,  data,  /No_Copy 

Handle_Free,  handle 
Return,  data 
End ;  of  picklist 

****************************************************************** 

9 

******************************************************************  « 


The  source  code  for  read._mos.pro  is  the  following: 

.  ***************************************************************** 

.  ***************************************************************** 

9 

Function  read_mos,  file,  imtype,  m,  n,  widl 

print, m,n 

;read  byte 

if  imtype  eq  0  then  begin 
openr,l,file 
newimg=bytarr(m,n) 
readu,l,newimg 
close,  1 

endif 
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;read  tiff 

if  imtype  eq  1  then  begin 

newimg=tiff_read(file) 

endif 

newimg=reverse(newimg,2);  reverses  the  rows 
return,  newimg 
end ;  of  read_mos 

9  ^ 

9 


The  source  code  for  seaminfo.pro  is  the  following: 

9  ****************^***********^i*************^:******^:*********:^***:^* 

9  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  H1  ^  ^  ^  ^  ^  ^  ^  ^  ^  H1  %  H1  ^  *£  ^  ^  ^  ^  ^  ^  % 

Pro  seaminfo_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET  JBUTTON’  THEN  BEGIN 
Widget_Control,  event.top,  /Destroy 

endif 

End ;  of  seaminfo_Event 

9  ^  ^  ^  ^  ^  H1  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

Function  seaminfo,  xoff,  yoff 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title-Define  Seam  Search’) 
filebase  =  Widget_Base(base,  Column=l) 

filelabel=  Widget_Label(base,  Value= ’Click  Mouse  At  Correlation  Line’) 
ok  =  Widget_Button(base,  Value=’OK’) 

Widget_Control,  base,  /Realize 
XManager,  ’seaminfo’,  base,  /Modal 
return, 0 

End ;  of  seaminfo 

9 

9  ^  ^  ^  ^  ^  ^  ^  ^ 
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APPENDIX  C 


Digital  Reconstruction  Program  and  Subprograms 

As  in  Appendix  B,  the  following  source  code  was  written  in  DDL.  The  source 
code  is  not  extensively  commented.  Its  inclusion  is  intended  mainly  for  reference  and 
documentation.  Some  of  the  source  code  was  obtained  while  at  a  training  class  given  by 
Research  Systems,  Inc.  During  programming,  it  was  decided  that  writing  one  large 
program  was  not  the  best  approach  because  there  would  be  a  desire  to  experiment  with 
variations  in  the  source  code.  Therefore,  different  functions  were  broken  out  into 
individual  programs  so  that  when  changes  were  made,  only  one  program  had  to  be 
modified.  The  programs  included  in  this  appendix  are  the  following: 

dhdouble.pro 

getd.pro 

getdim.pro 

getmisc.pro 

imcutarr.pro 

linelog.pro 

orient.pro 

read_dh.pro 

The  source  code  for  dhdouble.pro  is  the  following: 

. ***************************************************************** 

9 

.if:**************************************************************** 

pro  mask_event,  event 

Widget_Control,  event.top,  get_uvalue=info,  /No_Copy 
widget_controI,  info.holoid,  get_uvalue=holo,  /no_copy 
possibleEventTypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  NE  'DOWN'  THEN  RETURN 

info.mousex=event.x 

info.mousey=event.y 

;  Change  the  event  handler  for  the  draw  widget  and  turn  MOTION  events  ON. 
Widget_Control,  event. id,  Event_Pro='maskmove_event',  Draw_Motion_Events=l 

;  Take  over  color  index  1  for  the  zoom  box  drawing  color.  Store  the  current  (r,g,b)  values 
;  for  color  index  1  so  you  can  restore  the  current  colors  after  the  zoom  box  is  drawn. 
TVLct,  r,  g,  b,  /Get 
info.r_old  =  r(l) 
info.g_old  =  g(l) 
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info.b_old  =  b(l) 


widget_control,  info.holoid,  set_uvalue=holo,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end ;  of  mask_event 

.  ***************************************************************** 

pro  maskmove_event, event 

;  This  event  handler  continuously  draws  and  erases  the  image  box  until  it  receives  an  UP 
;  event  from  the  draw  widget.  Then  it  turns  draw  widget  motion  events  OFF  and  changes 
;  the  event  handler  for  the  draw  widget  back  to  zoom_PROCESS_EVENTS. 

;  Get  the  info  structure  out  of  the  top-level  base. 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 

widget_control,  info. draw  lid,  get_uvalue=recon,  /no_copy 

;  What  type  of  an  event  is  this? 

possibleEventT ypes  =  [  'DOWN',  'UP',  'MOTION',  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  EQ  'UP'  THEN  BEGIN 

;  If  this  is  an  UP  event,  you  need  to  erase  the  zoombox,  restore  the  user's  color 
;  table,  turn  motion  events  OFF,  set  the  draw  widget's  event  handler  back  to 
;  zoom_PROCESS_EVENTS,  and  draw  the  "zoomed"  plot  in  both  the  draw 
;  widget  and  the  pixmap. 

;  Turn  motion  events  and  button  events  off. 

Widget_Control,  event.id,  Draw_Motion_Events=0,  draw_button_events=0 

;  Restore  the  user's  color  table. 

TVLct,  info.r_old,  info.g_old,  info.b_old,  1 

info.maskxlow=info.mousex 

info.maskxhigh=event.x 

info.maskylow=info.mousey 

info.maskyhigh=event.y 

;  Make  sure  the  x  and  y  values  are  ordered  correctly 
IF  info.maskxlow  GT  info.maskxhigh  THEN  begin 
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info.maskxlow=event.x 

info.maskxhigh=info.mousex 

endif 

IF  info.maskylow  GT  info.maskyhigh  THEN  begin 
info.maskylow=event.y 
info.maskyhigh=info.mousey 

endif 

;  Scale  the  coordinates 

info.maskxlow=round(info.maskxlow*info.xscale) 

info.maskxhigh=round(info.maskxhigh*info.xscale) 

info.maskylow=round(info.maskylow*info.yscale) 

info.maskyhigh=round(info.maskyhigh*info.yscale) 

;  Display  the  masked  out  recon 

reeon(info.maskxlow:info.maskxhigh,info.maskylow:info.maskyhigh)=0 
if  info.disptype  eq  1  then  begin 
wset,info.widl 
erase 

tvscl,congrid(recon,5 1 2,5 1 2) 

wset,info.wid2 

erase 

tvscl,recon(info.zoomxlow:info.zoomxhigh,  $ 
info.zoomylow:info.zoomyhigh) 

endif 

if  info.disptype  eq  2  then  begin 
wset,info.widl 
erase 

tvscl,congrid(alog(recon  >  le-6),512,512) 

wset,info.wid2 

erase 

tvscl,alog(recon(info.zoomxlow:info.zoomxhigh,$ 
info.zoomylow:info.zoomyhigh)  >  le-6) 

endif 

widget_control,  info.drawlid,  set_uvalue=recon,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location  and  return. 

Widget_Control,  event. top,  Set_UValue=info,  /No_Copy 
RETURN 

ENDIF  ;  thisEvent  =  UP 

;  Most  of  the  action  in  this  event  handler  occurs  here  while  we  are  waiting  for  an  UP 
;  event  to  occur.  As  long  as  we  don't  get  it,  keep  erasing  the  old  zoom  box  and  drawing  a 
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;  new  one. 


;  Erase  the  old  zoom  box. 

WSet,  info.widl 

Device,  Copy  =  [0,  0,  512,  512,  0,  0,  info.boxwid] 

;  Update  the  dynamic  comer  of  the  zoom  box  to  the  current  cursor  location, 
newx  =  event.x 
newy  =  event.y 

;  Load  a  color  in  color  index  1  to  draw  the  zoom  box  with. 

TVLct,  255B,  OB,  OB,  1 

;  Draw  the  zoom  box. 

PlotS,  [info.mousex,  info.mousex,  newx,  newx,  info.mousex],  $ 

[info.mousey,  newy,  newy,  info.mousey,  info.mousey],  /Device,  Color=l 

widget_control,  info.drawlid,  set_uvalue=recon, /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end  ;  of  maskmove_event 

.  ******** ************************************************ ********* 


pro  zoom_event,  event 

Widget_Control,  event.top,  get_uvalue=info,  /No_Copy 
widget_control,  info.holoid,  get_uvalue=holo,  /no_copy 
possibleEventT ypes  =  [  DOWN’,  ’UP’,  ’MOTION’,  'SCROLL'  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  NE  ’DOWN’  THEN  RETURN 

info.mousex=event.x 

info.mousey=event.y 

;  Change  the  event  handler  for  the  draw  widget  and  turn  MOTION  events  ON. 
Widget_Control,  event.id,  Event_Pro='zoommove_event',  Draw_Motion_Events=l 

;  Take  over  color  index  1  for  the  zoom  box  drawing  color.  Store  the  current  (r,g,b)  values 

;  for  color  index  1  so  you  can  restore  the  current  colors  after  the  zoom  box  is  drawn. 

TVLct,  r,  g,  b,  /Get 

info.r_old  =  r(l) 

info.g_old  =  g(l) 

info.b_old  =  b(l) 
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widget_control,  info.holoid,  set_uvalue=holo,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end ;  of  zoom_event 

.  ***************************************************************** 


pro  zoommove_event, event 

;  This  event  handler  continuously  draws  and  erases  the  image  box  until  it  receives  an  UP 
;  event  from  the  draw  widget.  Then  it  turns  draw  widget  motion  events  OFF  and  changes 
;  the  event  handler  for  the  draw  widget  back  to  zoom_PROCESS_EVENTS. 

;  Get  the  info  structure  out  of  the  top-level  base. 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 

widget_control,  info.drawlid,  get_uvalue=recon,  /no_copy 

;  What  type  of  an  event  is  this? 

possibleEventT ypes  =  [  'DOWN',  ’UP’,  ’MOTION’,  ’SCROLL’  ] 
thisEvent  =  possibleEventTypes(event.type) 

IF  thisEvent  EQ  'UP'  THEN  BEGIN 

;  If  this  is  an  UP  event,  you  need  to  erase  the  zoombox,  restore  the  user’s  color 
;  table,  turn  motion  events  OFF,  set  the  draw  widget's  event  handler  back  to 
;  zoom_PROCESS_EVENTS ,  and  draw  the  "zoomed"  plot  in  both  the  draw 
;  widget  and  the  pixmap. 

;  Turn  motion  events  and  button  events  off. 

Widget_Control,  event.id,  Draw_Motion_Events=0,  draw_bntton_events=0 

;  Restore  the  user's  color  table. 

TVLct,  info.r_old,  info.g_old,  info.b_old,  1 

info.zoomxlow=info.mousex 

info.zoomxhigh=event.x 

info.zoomylow=info.mousey 

info.zoomyhigh=event.y 

;  Make  sure  the  x  and  y  values  are  ordered  correctly. 

IF  info.zoomxlow  GT  info.zoomxhigh  THEN  begin 
info.zoomxlow=event.x 
info.zoomxhigh=info.mousex 
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endif 

IF  info.zoomylow  GT  info.zoomyhigh  THEN  begin 
info.zoomylow=event.y 
info.zoomyhigh=info.mousey 

endif 

;  Scale  the  coordinates 

info.zoomxlow=round(info.zoomxlow*info.xscale) 

info.zoomxhigh=round(info.zoomxhigh*info.xscale) 

info.zoomylow=round(info.zoomylow*info.yscale) 

info.zoomyhigh=round(info.zoomyhigh*info.yscale) 

print,info 

;  Display  the  zoomed  recon 

wset,  info.widl 

erase 

wset,  info.wid2 
erase 

if  info.disptype  eq  1  then  begin 
wset,info.widl 
tvscl,congrid(recon,512,512) 
wset,info.wid2 

tvscl,congrid(recon(info.zoomxlow:info.zoomxhigh,$ 
info.zoomylow:info.zoomyhigh),  1 50, 1 50) 

endif 

if  info.disptype  eq  2  then  begin 
wset, info.widl 

tvscl,congrid(alog(recon  >  le-6),512,512) 
wset,info.wid2 

tvscl,congrid(alog(recon(info.zoomxlow:info.zoomxhigh,$ 
info.zoomylow:info.zoomyhigh)  >  le-6),150,150) 

endif 

widget_control,  info,  draw  lid,  set_uvalue=recon,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location  and  return. 
widget_Control,  event.top,  Set_UValue=info,  /No_Copy 
RETURN 

ENDIF  ;  thisEvent  =  UP 

;  MoSt  of  the  action  in  this  event  handler  occurs  here  while  we  are  waiting  for  an  UP 
;  event  to  occur.  As  long  as  we  don't  get  it,  keep  erasing  the  old  zoom  box  and  drawing  a 
;  new  one. 
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;  Erase  the  old  zoom  box. 

WSet,  info.widl 

Device,  Copy  =  [0,  0,  512,  512,  0, 0,  info.boxwid] 

;  Update  the  dynamic  comer  of  the  zoom  box  to  the  current  cursor  location, 
newx  =  event,  x 
newy  =  event.y 

;  Load  a  color  in  color  index  1  to  draw  the  zoom  box  with. 

TVLct,  255B,  OB,  OB,  1 

;  Draw  the  zoom  box. 

PlotS,  [info.mousex,  info.mousex,  newx,  newx,  info.mousex],  $ 

[info.mousey,  newy,  newy,  info.mousey,  info.mousey],  /Device,  Color=l 

widget_control,  info.drawlid,  set_uvalue=recon,  /no_copy 

;  Put  the  info  structure  back  into  its  storage  location. 

Widget_Control,  event.top,  Set_UValue=info,  /No_Copy 

end ;  of  zoommove_event 

.  ***************************************************************** 


PRO  dhdouble_mask_Event,  event 

widget_control, event.top, get_uvalue=info 
print,  info 

widget_control,event.id,get_uvalue=value 
widget_control,info.holoid,get_uvalue=holo 
help,  /memory 

CASE  value  OF 

'DIM'  :  begin  ;obtain  image  X  &  Y  size 

widget_control, event.top, tlb_get_offset=offsets 
imdim=getdim(offsets(0)+50,offsets(l)+75,  info.m,  info.n) 
if  imdim.xsize  ne  0  then  begin 
info.m=imdim.xsize 
info.n=imdim.ysize 
info.mc=imdim.xsize 
info.nc=imdim.ysize 
info.xscale=double(info.mc/512.) 
info.yscale=double(info.nc/5 1 2.) 
print, 'info.yscale',info.yscale,size(info.yscale) 
widget_control,event.top,set_uvalue=info 

endif 
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end 

'ORIENT' :  begin  ; obtain  orientation  of  plane  wave 

widget_control,event.top,tlb_get_offset=offsets 

orientdim=orient(offsets(0)+50,offsets(l)+75,info.theta,info.phi) 

info.theta=double(orientdim.theta) 

info.phi=double(orientdim.phi) 

widget_control,event.top,set_uvalue=info 

end 

'D'  :  begin  ;obtain  distance  between  holo  and  object 
widget_control,event.top,tlb_get_offset=offsets 
info.d=getd(offsets(0)+50,offsets(l)+75,info.d) 
info.d=double(info.d) 
widget_control,event.top,set_uvalue=info 
end 

'MISC'  :  begin 

widget_controI,event.top,tlb_get_offset=offsets 
misc=getmisc(offsets(0)+50,offsets(l)+75,info.dx,info.dy  $ 
,info.lambda,info.mag) 
info.dx=double(misc.dx) 
info.dy=double(misc.dy) 
info.lambda=double(misc.lambda) 
info.mag=double(misc.mag) 
widget_control,event.top,set_uvalue=info 
end 

'LOAD'  :begin  ;read  interference  image 

widget_control,event.top,tlb_get_offset=offsets 
pick=dialog_pickfile(/read,  path=info.impath,  get_path=impath) 
info.impath=impath 
info.pick=pick 

;determine  the  image  type 

pick=strlowcase(pick) 

type=strmid(pick,strlen(pick)-3,3) 

if  type  eq  'byt'  then  imtype=0 
if  type  eq  'tif  then  imtype=l 
if  type  eq  'byt'  or  type  eq  'tif  then  begin 
info .  imtype=imtype 

holo=read_dh(info.pick,info.imtype,  info.m,  info.n,  info.widl) 
widget_control,  event.top,  set_uvalue=info 
widget_control,  info.holoid,  set_uvalue=holo,  /no_copy 

endif 

end 

'RELOAD'  : begin 
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;  the  read  program  reads  the  current  file  and  displays  it,  it  returns  holo 
holo=read_dh(info.pick,info.imtype,  info.m,  info.n,  info.widl) 

info.mc=info.m 

info.nc=info.n 

widget_control,event.top,set_uvalue=info 
widget_control,  info.holoid,  set_uvalue=holo,  /no_copy 
end 

'CUT'  :begin 

widget_control,event.top,tlb_get_offset=offsets 

;  send  me  &  nc  so  that  present  values  can  be  displayed 
cutdata=imcutarr(offsets(0)+50,  offsets(l)+75,  info.mc,  info.nc,  $ 
info.pick,info.imtype,  info.m,  info.n,  info.widl) 

if  cutdata.accept  eq  0  then  begin  ;use  current  dimensions 
info.mc=info.mc 
info.nc=info.nc 

endif 

if  cutdata.accept  eq  1  then  begin;  if  not  1,  then  rest  is  skipped 
info.mc=cutdata.newxsize 
info.nc=cutdata.newysize 
info.xscale=double(info.mc/512.) 
info.yscale=double(info.nc/512.) 
if  cutdata.newcenter  eq  1  then  begin 

cursor,  subx,  suby,  /device,  /down 

info.suby=suby 

info.subx=subx 

endif 

holo=extrac(holo,info.subx*info.m/5 12-info.mc/2,  $ 
info.suby*info.n/512-info.nc/2,  info.mc,  info.nc) 
widget_control,info.holoid,set_uvalue=holo 
wset,  info.widl 
erase 

tvscl,congrid(holo,5 1 2,5 1 2) 

endif 

widget_control,event.top,set_uvalue=info 

end 

’RELOADCUT'  :begin 
;read  byte 

if  info.imtype  eq  0  then  begin 
openr,l,info.pick 
holo=bytarr(info.m,  info.n) 
readu,l,holo 


close,  1 


endif 
;read  tiff 

if  info.imtype  eq  1  then  begin 

holo=tiff_read(info.pick) 

endif 

;  extract  the  subimage  from  the  original  image 

holo=extrac(holo,info.subx*info.m/512-info.mc/2,  $  \ 

info.suby*info.n/512-info.nc/2,info.mc,info.nc) 

wset,  info.widl 
erase 

tvscl,congrid(holo,5 1 2,5 1 2) 

widget_control,  info.holoid,  set_uvalue=holo,  /no_copy 
end 

'SAVE'  : begin 

widget_control,info.draw  1  id,get_uvalue=recon 
if  info.disptype  eq  2  then  begin 

recon=alog(double(recon)+ 1  d-6) 

endif 

recon=bytscl(recon) 

pick=dialog_pickfile(/write,  path=info.impath,  get_path=impath) 

info.impath=impath 

tiff_write,pick, recon 

widget_control,event.top,set_uvalue=info 

end 

’MASK'  : begin 

widget_control,  info.drawlid,  draw_Button_Events=l,  $ 
event_pro='mask_event’ 
wset,info.boxwid 

device,  copy  =  [0, 0,  512,  512,  0, 0,  info.widl] 
widget_control,  info.drawlid,  set_uvalue=recon,  /no_copy 
widget_Control,  event.top,  set_uvalue=info 
end 

'LINELOG' :  begin  > 

widget_control, event.top, tlb_get_offset=offsets 
widget_control,info.draw  1  id,get_uvalue=recon 
disptype=linelog(offsets(0)+50,  offsets(  1  )+75) 
if  disptype  eq  1  then  begin 
wset, info.widl 
erase 

tvscl,congrid(recon,5 12,512) 
wset,info.wid2 
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erase 

tvscl,recon(info.zoomxlow:info.zoomxhigh,  $ 
info.zoomylow:info.zoomyhigh) 

endif 

if  disptype  eq  2  then  begin 
wset,info.widl 
erase 

tvscl,congrid(alog(double(recon)+ld-6),512,5 12) 

wset,info.wid2 

erase 

tvscl,alog(double(recon(info.zoomxlow:info.zoomxhigh,  $ 
info.zoomylow:info.zoomyhigh))+ 1  d-6) 

endif 

info.disptype=disptype 

widget_control,event.top,set_uvalue=info 

end 

'RECON' : begin 

calculate  planewave  and  gaussian 
k=2.*  !dpi/info.lambda 

kx=k*cos(info. theta*  !dpi/l  80d)*sin(info. phi*  !dpi/180d) 

ky=k*sin(info.theta*!dpi/180d)*sin(info.phi*!dpi/180d) 

vectx=(dindgen(info.mc)-round(info.mc/2))*info.dx/info.mag 

arx=transpose(make_array(info.nc,  1  ,/double,value=  1 .0)) 

vecty=(transpose(dindgen(info .  nc)-round(info  .nc/2)))  *info.dy/info .  mag 

ary=make_array(info.mc,  1  ,/double,value=  1 .0) 

planewave=kx  *  (vectx#arx)+ky*(ary#vecty) 

planewave=exp(dcomplex(0,temporary(planewave))) 

gauss=-((vectx#arx)A2+(ary#vecty)A2)*k/2.d/info.d 

gauss=exp(dcomplex(0,temporary(gauss))) 

recon=double(holo)*gauss*planewave 

gauss=0 

planewave=0 

;  Positive  #  indicates  an  inverse  Fourier  transform 
recon=fft(temporary(recon),- 1  ,double=l ,  /overwrite) 

recon=abs(temporary(recon))A2 

recon=shift(temporary(recon),info.mc/2,info.nc/2) 

beep 

if  info.disptype  eq  1  then  begin 
wset,info.widl 
tvscl,congrid(recon,5 12,5 12) 
wset,info.wid2 
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tvscl,congrid(recon(info.zoomxlow:info.zoomxhigh,  $ 
info.zoomylow.info.zoomyhigh),  1 50, 1 50) 

endif 

if  info.disptype  eq  2  then  begin 
wset,info.widl 

tvscl,congrid(alog(double(recon)+ 1  d-6),5 12,512) 
wset,info.wid2 

tvscl,congrid(alog(double(recon(info.zoomxlow:info.zoomxhigh,  $ 
info.zoomylow:info.zoomyhigh))+ld-6),150,150) 

endif 

widget_control,event.top,set_uvalue=info 
widget_control,info.holoid,set_uvalue=holo,  /no_copy 
widget_control,info.draw  1  id,set_uvalue=recon,  /no_copy 
end 

'ZOOM'  rbegin 

widget_control,  info.drawlid,  draw_Button_Events=l,  $ 
event_pro=’zoom_event' 
wset,info.boxwid 

device,  copy  =  [0,  0,  512,  512,  0,  0,  info.widl] 
widget_control,  info.drawlid,  set_uvalue=recon,  /no_copy 
widget_Control,  event.top,  set_uvalue=info 
end 

'EXIT'  :  Widget_Control,  event.top,  /Destroy 
ENDCASE 

END  ;  of  dhdouble_mask_Event 

•  5jf  ^  *1^  ^  ^  ^  ^  ^  ^  ^J>  ^  vj^  ^  ^  v|#  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  vj>  %|«  »|»  ^ 

PRO  dhdouble 

;  images  are  transferred  from  bottom  to  top,  i.e.  the  row  with  a  0  subscript  is  drawn  on 

;  the  bottom 

!order=0 

device,  retain=2 

base  =  Widget_Base(Title='Digital  Holography  (double  precision)',  Column=2) 
drawbase=widget_base(base,column=2) 
btnbase=widget_base(base,column=  1 ) 

btn  1  =widget_button(btnbase,value=Tnput  Image  Dimensions’, uvalue='DIM') 
btn2=widget_button(btnbase,value='Input  Plane  Orientation', uvalue='ORIENT’) 
btn3=widget_button(btnbase,value='Input  D',uvalue='D') 
btnl3=widget_button(btnbase,value='Input  dx,  dy,  lambda,  mag',uvalue=’MISC’) 
btn4=widget_button(btnbase,value='Load  Interference  Image', uvalue='LO AD') 
btn5=widget_button(btnbase,value='Re-Load  Interference  Image', uvalue='RELOAD') 
btn6=widget_button(btnbase,value= 'Extract  Partial  Image',uvalue='CUT') 
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btnl  l=widget_button(btnbase,value='Re-Load  Extracted  Image',uvalue='RELOADCUT) 
btn7=widget_button(btnbase,value='Reconstruct  Hologram' ,uvalue= 'RECON') 
btn  1 2=widget_button(btnbase,value='Linear/Log  Display' ,uvalue='LINELOG') 
btnl4=widget_button(btnbase,value='Zoom’,uvalue='ZOOM') 
btnlO=widget_button(btnbase,value='Mask  Image  Portion’, uvalue='MASK’) 
btn8=widget_button(btnbase,value='Save  Reconstructed  Image  (TIF),,uvalue='SAVE') 
btn9=widget_button(btnbase,value=’Exit',uvalue='EXIT') 
draw  1  =widget_draw(drawbase,xsize=5 1 2,ysize=5 1 2) 
draw2=widget_draw(btnbase,xsize=  1 50,ysize=  1 50) 

Widget_Control,  base,  /Realize 

widget_control,draw  1  ,get_value=wid  1 

widget_control,draw2,get_value=wid2 

Window,  /Free,  XSize=512,  YSize=512,  /Pixmap 

boxwid  =  !D  Window 

xscale=  1 024d/5 1 2d 

yscale= 1 024d/5 1 2d 

info={m:1024,  n:1024,  mc:1024,  nc:1024,  theta:0d,phi:0d,  d:0.218d,  $ 

widl:widl,  wid2:wid2,  holoid:drawbase,  draw  lid:  drawl,  imtype:0,  $ 
disptype:2,  subx:1281,  suby:1281,  pick:"",  boxwid:boxwid,  mousex:0,  $ 
mousey:0,  r_old:0,  g_old:0,  b_old:0,  maskxlow:0,  maskxhigh:0,  $ 
maskylow:0,  maskyhigh:0,  xscale:xscale,  yscale:yscale,  impath:'e:Y,  $ 
mag:7.4664d,  lambda:632.8d-9,  dx:6.8d-6,  dy:6.8d-6,  zoomxlow:0,  $ 
zoomxhigh:127,  zoomylow:0,  zoomyhigh:127} 
holo=0b 
recon=0b 
loadct,0 
wset,widl 
erase 

widget_control,base,set_uvalue=info 
widget_control,drawbase,set_uvalue=holo 
widget_control,draw  1  ,set_uvalue=recon 

XManager,  'dhdouble_mask',  base,  Event='dhdouble_mask_Event' 

END  ;  of  PRO  dhdouble 

.***************************************************************** 

? 

.  ***************************************************************** 


The  source  code  for  getd.pro  is  the  following: 

.  ***************************************************************** 
****************************************************************** 


Pro  getd_Event,  event 
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name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  rWIDGET_BUTTON'  THEN  BEGIN 

Widget_ControI,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

'CANCEL'  :  Begin 

data  =  info.currentd 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

'ACCEPT'  :  Begin 

Widget_Control,  info.d_id,  Get_Value=d 
data  =  d 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

ENDCASE 


ENDIF 

END  ;  of  getd_Event 
Function  getd,  xoff,  yoff,  currentd 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title=’Enter  D  Value') 

filebase  =  Widget_Base(base,  Column=l) 

filelabel  =  Widget_Label(filebase,  Value=’Input  D') 

d_id  =  CW_Field(base,  Title='D: ',  Value=currentd,  /float,  /Retum_Events) 

btnbase  =  Widget_Base(base,  Row=l,  /Frame) 

cancel  =  Widget_Button(btnbase,  Value='CANCEL') 

done  =  Widget_Button(btnbase,  Value='  ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {d_id:d_id,  handle:handle,currentd:currentd} 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'getd',  base,  /Modal 

;Must  use  handle  because  all  widgets  are  destroyed  after  Xmanager. 
Handle_Value,  handle,  data,  /No_Copy 


Handle_Free,  handle 
Return,  data 
End  ;  Function  getd 
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.  ****************************************** *********************** 
> 

.  ***************************************************************** 


The  source  code  for  getdim.pro  is  the  following: 

.  ***************************************************************** 
> 

.  ***************************************************************** 
5 

Pro  getdim_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  ’WIDGET_BUTT ON'  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

'CANCEL'  :  Begin 

data={xsize:0} 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

'ACCEPT'  :  Begin 

Widget_Control,  info.xsize_id,  Get_Value=xsize 
Widget_Control,  info.ysize_id,  Get_Value=ysize 
data=  {xsizerxsize,  ysize:ysize} 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

ENDCASE 


ENDIF 

END  ;  of  getdim_Event 

.  ************* **************************************************** 


Function  getdim,  xoff,  yoff,  currentx,  currenty 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  Title= 'Enter  File  Sizes') 
filebase  =  Widget_Base(base,  Column=l) 

filelabel  =  WidgetJLabel (filebase,  Value='Input  Image  Dimensions') 

xsizejd  =  CW_Field(base,  Title='X  Size:  ’,  Value=currentx,  /Integer,  /Retum_Events) 

ysize_id  =  CW_Field(base,  Title='Y  Size: ',  Value=currenty,  /Integer,  /Retum_Events) 

btnbase  =  Widget JBase(base,  Row=l,  /Frame) 

cancel  =  Widget_Button(btnbase,  Value='CANCEL') 
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done  =  Widget_Button(btnbase,  Value-  ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {xsize_id:xsize_id,  ysize_id:ysize_id,  handlerhandle,  currentxrcurrentx,  $ 
currenty:currenty } 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'getdim',  base,  /Modal 
Handle_Value,  handle,  data,  /No_Copy 
Handle_Free,  handle 
Return,  data 

End  ;  of  Function  getdim 

.  ***************************************************************** 
. ***************************************************************** 


The  source  code  for  getmisc.pro  is  the  following: 

.  sjcjjc^ilcjJjjfc*********************************************************** 

.  ^c^c^c^c^c^c^^c^c^c^c^c^c^c^:^:^:^:^^:^:^:^****************************************** 
) 

Pro  getmisc_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET_BUTTON’  THEN  BEGIN 

Widget_Control,  event. top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event. id,  Get_Value=test 

CASE  test  OF 

'CANCEL'  :  Begin 

data=  { dx:currentdx,dy:currentdy,lambda:currentlambda,  $ 
magrcurrentmag } 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

'ACCEPT'  :  Begin 

Widget_Control,  info.dx_id,  Get_Value=dx 
Widget_Control,  info.dy_id,  Get_Value=dy 
Widget_Control,  info.lambda_id,  Get_Value=lambda 
Widget_Control,  info.mag_id,  Get_Value=mag 
data  =  {dx:dx,dy:dy,lambda:lambda,mag:mag} 
Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
End 

ENDCASE 
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ENDEF 

END  ;  of  getmisc_Event 

.  ***************************************************************** 


Function  getmisc,  xoff,  yoff,  currentdx,  currentdy,  currentlambda,  currentmag 
base  =  Widget_Base(Column=  1 ,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Miscellaneous  Parameters') 
filebase  =  Widget_Base(base,  Column=l) 

{  filelabel  =  Widget_Label(filebase,  Value='Input  Parameters’) 

dx_id  =  CW_Field(base,  Title='dx: Value=currentdx,  /float,  /Retum_Events) 
dy_id  =  CW_Field(base,  Title='dy: Value=currentdy,  /float,  /Retum_Events) 
r  lambdajd  =  CW_Field(base,  Title='Lambda: ',  Value=currentlambda,  $ 

/float,  /Retum_Events) 

mag_id  =  CW_Field(base,  Title=’Mag:  ’,  Value=currentmag,  /float,  /Retum_Events) 
btnbase  =  Widget_Base(base,  Row=  1 ,  /Frame) 
cancel  =  Widget_Button(btnbase,  Value='CANCEL') 
done  =  Widget_Button(btnbase,  Value=' ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {dx_id:dx_id,  dy_id:dy_id,  lambda_id:lambda_id,  mag_id:mag_id,  handle:handle) 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 

XManager,  'getmisc',  base,  /Modal 

Handle_Value,  handle,  data,  /No_Copy 

Handle_Free,  handle 

Return,  data 

End  ;  of  Function  getmisc 

.  ***************************************************************** 

. ***************************************************************** 

5 

The  source  code  for  imcutarr.pro  is  the  following: 

.***************************************************************** 

.if;**************************************************************** 

i 

Pro  imcutarr_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET_BUTTON'  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 
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'ACCEPT'  :  Begin 

Widget_Control,  info.newxsize_id,  Get_Value=newxsize 
Widget_Control,  info.newysize_id,  Get_Value=newysize 
data  =  {newxsize:newxsize,newysize:newysize,  $ 
accept:  1  ,newcenter:  1 } 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 

widget_Control,  info.base2_id,  /realize 
Widget_Control,  info.base2_id,  Set_UValue=info 
XManager,  'redisplay',  info.base2_id,  /Modal 
End 

'CANCEL'  :  Begin 
data  =  {accept:0} 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 
End 

ENDCASE 


ENDIF 

END  ;  of  imcutarr_Event 

^  ^  ^  ^  v|#  ^  ^  ^  ^  ^  v|/  ^  ^  ^  sj#  ^  ^  ^  jJ/  vj/  ^  jJj  ^ 

Pro  redisplay_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  ’WIDGET.BUTTON’  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event. id,  Get_Value=test 

CASE  test  OF 

'NO'  :  Begin 

Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 

widget_Control,  info.base3_id,  /realize 
Widget_Control,  info.base3_id,  Set_UValue=info 
XManager,  'mouseinfo',  info.base3_id,  /Modal 
end 

'YES'  :  Begin 

holo=read_dh(info.file,  info.imtype,  info.m,  info.n,  info.widl) 
Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 

widget_Control,  info.base3_id,  /realize 
Widget_Control,  info.base3_id,  Set_UValue=info 
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XManager,  'mouseinfo',  info.base3_id,  /Modal 
end 

ENDCASE 


endif 

End  ;  redisplay_Event 

.  ************************** **************************** *********** 

Pro  mouseinfo_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  ’WIDGET.BUTTON’  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Handle_Value,  info.handle,  data,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

’CANCEL’  :  Begin 
data.newcenter=0 

Handle_VaIue,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 
end 

'OK'  :  Begin 

data.newcenter=l 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
wshow,  info.widl 
end 

ENDCASE 


endif 

End  ;  mouseinfo_Event 

.  ***************************************************************** 

Function  imcutarr,  xoff,  yoff,  currentxsize,  currentysize,  file,  imtype,  m,  n,  widl 

base  =  Widget_Base(Column=  1 ,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Cut  Interference  Image') 
filebase  =  Widget_Base(base,  Column=l) 

filelabel  =  Widget_Label(filebase,  Value='Specify  New  Image  Dimensions') 
base3  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title=’Select  Image  Center’) 
filebase3  =  Widget_Base(base3,  Column=l) 

filelabel3=  Widget_Label(base3,  Value='Click  Mouse  At  Desired  Center') 
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ok  =  Widget_Button(base3,  Value='OK') 
can=Widget_Button(base3,  Value=’CANCEL') 

base2=Widget_Base(Column=  1 ,  XOffset=xoff,  YOffset=yoff,  Title='Redisplay') 
filelabel2=Widget_Label(base2,  Value='Redisplay  Original  Image?') 
yes=Widget_Button(base2,  Value='YES') 
no=W idget_Button(base2,  V alue='N O’) 

newxsizejd  =  CW_Field(base,  Title='New  X  Size: Value=currentxsize,  $ 

/long,  /Retum_Events) 

newysize_id  =  CW_Field(base,  Title=’New  Y  Size: Yalue=currentysize,  $ 

/long,  /Retum_Events) 

btnbase  =  Widget_Base(base,  Row=l,  /Frame) 
cancel  =  Widget_Button(btnbase,  Value='CANCEL') 
done  =  Widget_Button(btnbase,  Value=' ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {newxsize_id:newxsize_id,newysize_id:newysize_id,handle:handle,  $ 
base2_id:base2,  base3_id:base3,currentxsize:currentxsize,  $ 
currentysizexurrentysize,  flle:flle,  imtype:imtype,  m:m,  n:n,  widl:widl } 
Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'imcutarr',  base,  /Modal 
Handle_Value,  handle,  data,  /No_Copy 
Handle_Free,  handle 
Return,  data 

End  ;  of  Function  imcutarr 

.  ***************************************************************** 
? 

.***************************************************************** 


The  source  code  for  linelog.pro  is  the  following: 

.******** ********* ****** ********************* ********************* 
9 

.  ***************************************************************** 
9 

Pro  linelog_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 

IF  name  EQ  'WIDGET_BUTTON'  THEN  BEGIN 

Widget_Control,  event,  top,  Get_UValue=handle 
Widget_Control,  event. id,  Get_Value=test 

CASE  test  OF 

'LINEAR'  :Begin 

Handle_Value,  handle,  1,  /Set 
Widget_Control,  event.top,  /Destroy 


End 

'LOG'  :Begin 

HandleJValue,  handle,  2,  /Set 
Widget_Control,  event.top,  /Destroy 
End 

'CANCEL'  :Begin 

Handle_Value,  handle,  3,  /Set 
Widget_Control,  event.top,  /Destroy 
End 

ENDCASE 


ENDEF 

END  ;  of  Pro  linelog 

.  ***************************************************************** 
? 

function  linelog,  xoff,  yoff 

base  =  WidgetJBase(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Linear/Log  Display') 
filebase  =  Widget_Base(base,  Column=l) 

filelabel  =  Widget_Label(filebase,  Value='Select  Display  Method') 
linear  =  Widget_Button(base,  Value-LINEAR') 
log  =  Widget_Button(base,  Value='LOG') 
cancel  =  Widget_Button(base,  Value='CANCEL') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

Widget_Control,  base,  Set_UValue=handle 

XManager,  'linelog',  base,  /Modal 

Handle_Value,  handle,  data,  /No_Copy 

Handle_Free,  handle 

Return,  data 

End  ;  function  linelog 

.  ************ ******************************** ********************* 

> 

.  ***************************************************************** 


The  source  code  for  orient.pro  is  the  following: 

.  ***************************************************************** 

.  ^**** ******************** **************************************** 
y 

Pro  orient_Event,  event 

name  =  Tag_Names(event,  /Structure_Name) 


IF  name  EQ  ’WIDGET_BUTTON'  THEN  BEGIN 

Widget_Control,  event.top,  Get_UValue=info,  /No_Copy 
Widget_Control,  event.id,  Get_Value=test 

CASE  test  OF 

’CANCEL'  :  Begin 

data  =  {theta:currenttheta,phi:currentphi} 
HandleJValue,  info.handle,  data,  /Set,  /No_Copy 
Widget_Control,  event.top,  /Destroy 
End 

’ACCEPT’  :  Begin 

Widget_Control,  info.theta_id,  Get_Value=theta 
Widget_Control,  info.phi_id,  Get_Value=phi 
data=  {theta:theta,phi:phi} 

Handle_Value,  info.handle,  data,  /Set,  /No_Copy 

Widget_Control,  event.top,  /Destroy 

End 

ENDCASE 


ENDIF 

END  ;  of  orient_Event 

Function  orient,  xoff,  yoff,  currenttheta,  currentphi 

base  =  Widget_Base(Column=l,  XOffset=xoff,  YOffset=yoff,  $ 

Title='Plane  Wave  Orientation') 
filebase  =  Widget_Base(base,  Column=l) 

filelabel  =  Widget_Label(filebase,  Value='Input  Plane  Orientation  (degrees)') 
theta_id  =  CW_Field(base,  Title='Theta: ',  Value=currenttheta,  /float,  /Retum_Events) 
phi_id  =  CW_Field(base,  Title='Phi: ',  Value=currentphi,  /float,  /Retum_Events) 
btnbase  =  Widget_Base(base,  Row=l,  /Frame) 
cancel  =  Widget_Button(btnbase,  Value='CANCEL') 
done  =  Widget_Button(btnbase,  Value=' ACCEPT') 

Widget_Control,  base,  /Realize 
handle  =  handle_create() 

info  =  {theta_id:theta_id,  phi_id:phi_id,  handle:handle} 

Widget_Control,  base,  Set_UValue=info,  /No_Copy 
XManager,  'orient',  base,  /Modal 
Handle_Value,  handle,  data,  /No_Copy 
Handle_Free,  handle 
Return,  data 

End  ;  of  Function  orient 


.  ***************************************************************** 


The  source  code  for  read_dh.pro  is  the  following: 

.***************************************************************** 

9 

.***************************************************************** 

9 

Function  read_dh,  file,  imtype,  m,  n,  widl 
;read  byte 

if  imtype  eq  0  then  begin 
openr,l,file 
holo=bytarr(m,n) 
readu,l,holo 
close,  1 

endif 
;read  tiff 

if  imtype  eq  1  then  begin 
holo=tiff_read(file) 
print, size(holo) 

endif 

holo=reverse(holo,2);  reverses  the  rows 

wset,  widl 

erase 

tvscl,  congrid(holo,512,512) 
return,  holo 

end  ;  of  Function  read_dh 

.  ***************************************************************** 
9 

. ***************************************************************** 


( 
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APPENDIX  D 


Relative  Error  Propagation 


A  key  element  of  interpreting  the  experimental  data  is  error  analysis.  This  section 
will  be  devoted  to  a  short  summary  of  error  analysis,  and  in  particular,  relative  error 
propagation.  First  of  all,  relative  error,  £%. ,  is  defined  as  [Stoer,  Bulirsch,  1993,  p.  12] 


£7,  = 


Xi  Xj 


(D.l) 


where  x)  is  the  estimate  or  measured  value  of  Xj ,  and  Xj  is  the  true  value.  Once  the 
relative  errors  of  all  variables  are  known,  the  total  relative  error  must  be  determined.  This 
total  relative  error  depends  on  the  mathematical  operations  performed  in  reaching  the 
final  result. 


k 


The  equation  which  applies  in  determining  relative  error  propagation  resulting 
from  a  mathematical  operation  is  given  by  the  first  order  approximation  [Stoer,  Bulirsch, 
1993,  p. 13] 


y  xj  00  c 

dxj  Xl 


(D.l) 


For  example,  if 


<t>{x,y)  —  x-y ,  (D.3) 

where  the  variable  x  has  a  relative  error  of  £x  and  the  variable  y  has  a  relative  error  of 
ey ,  then  the  relative  of  (j)(x,y )  is  given  by  [Stoer,  Bulirsch,  1993,  p.  14] 


£<p(x,y) 


x  d(xy) 

- 3  e 

xy  ox 

1  1 

-y£x  +~ 


x  + 

X£y 


y  d(xy) 

xy  dy 


(D.4) 


-  £x  +  £y- 

In  a  similar  fashion,  the  following  relative  error  propagation  formulas  for  common 
mathematical  operations  can  be  computed  from  application  of  Eq.  (D.2)  [Stoer,  Bulirsch, 
1993,  p.  14]: 


(D.5) 

<j>(x,y)  =  x+y=*£^Xty)=  £x+  y  £y 

x  l  y  x  1  y 

(D.6) 

<P(x, y)  =  x  y=*£(/>(X>y)  =  x_y£x  xy_y£y 

(D.7) 

Care  must  be  taken  in  applying  Eq.  (D.5).  At  first,  it  may  seem  that  the  relative  errors  are 
simply  subtracted.  However,  this  equation  actually  is 


r 
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One  final  point  must  be  made  regarding  relative  errors.  The  point  is  that  in  order 
to  calculate  the  relative  error  as  defined  in  Eq.  (D.l),  the  true  value  of  the  variables  must 
be  known.  Obviously,  if  the  true  values  are  known,  then  there  are  no  errors  to  be 
concerned  about.  Therefore,  a  method  of  estimating  the  relative  error  must  be  used.  A 
reasonable  approximation  to  use  is  to  use  the  estimate  or  measured  value  of  the  variables 
in  the  denominator  of  Eq.  (D.l).  In  addition,  the  numerator  of  Eq.  (D.l),  known  as  the 
absolute  error,  is  found  by  applying  a  reasonable  estimate  of  how  close  the  measured 
value  is  to  the  true  value.  A  good  example  is  found  in  the  case  of  measuring  a  length 
with  a  ruler.  Generally,  it  is  reasonable  to  assume  that  a  measurement  can  be  made  to 
within  the  length  of  the  smallest  scale  on  the  ruler.  For  example,  if  the  smallest  scale  is 
1mm,  then  a  good  estimate  is  that  the  true  value  should  be  within  ±  0.5mm. 
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